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ABSTRACT 

We have carried out the first survey of the pulsational line profile variability in rapidly os- 
cillating Ap (roAp) stars. We analysed high signal-to-noise time-series observations of ten 
sharp-lined roAp stars obtained with the high-resolution spectrographs attached to the VLT 
and CFHT telescopes. We investigated in detail the variations of Pr III, Nd II, Nd III and Tb III 
lines and discovered a prominent change of the profile variability pattern with height in the 
atmospheres of all studied roAp stars. We show that, in every investigated star, profile vari- 
ability of at least one rare-earth ion is characterized by unusual blue-to-red moving features, 
which we previously discovered in the time-resolved spectra of the roAp star 7 Equ. This 
behaviour is common in rapidly rotating non-radial pulsators but is inexplicable in the frame- 
work of the standard oblique pulsator model of slowly rotating roAp stars. Using analysis of 
the line profile moments and spectrum synthesis calculations, we demonstrate that unusual os- 
cillations in spectral lines of roAp stars arise from the pulsational modulation of line widths. 
This variation occurs approximately in quadrature with the radial velocity changes, and its 
amplitude rapidly increases with height in stellar atmosphere. We propose that the line width 
modulation is a consequence of the periodic expansion and compression of turbulent layers 
in the upper atmospheres of roAp stars. Thus, the line profile changes observed in slowly ro- 
tating magnetic pulsators should be interpreted as a superposition of two types of variability: 
the usual time-dependent velocity field due to an oblique low-order pulsation mode and an 
additional line width modulation, synchronized with the changes of stellar radius. Our expla- 
nation of the line profile variations of roAp stars solves the long-standing observational puzzle 
and opens new possibilities for constraining geometric and physical properties of the stellar 
magnetoacoustic pulsations. 
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1 INTRODUCTION 

About 10 % to 20 % of upper main sequence stars are characterized 
by remarkably rich line spectra, often containing numerous uniden- 
tified features. Compared to the solar case, overabundances of up 
to a few dex are often inferred for some iron-peak and rare-earth el- 
ements, wher eas some other chemical elements are found to be un- 
derabundant dRvabchikova e t al. 2004). Some of these chemically 
peculiar (Ap) stars also exhibit organized magnetic fields with a 
typical strength of a few kG. Chemical peculiarities are believed to 

* Based on observations collected at the European Southern Observatory 
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trieved from the ESO Science Archive, 
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result from the influence of the magnetic field on the diffusing ions, 
possibly in combin ation with the influence of a weak, magnetically- 
directed wind (e.g.. lBabellll992h . 

More than 30 cool Ap stars exhibit high-overtone, low-degree, 
non-radia l p-mode pulsations wi th periods in the range of 6-21 
minutes (Kurtz & Martinez 2000). These rapidly oscillating Ap 
(roAp) stars are key objects for asteroseismology, which presently 
is the most powerful tool for testing theories of stellar structure 
and evolution. The observed pulsation amplitudes of roAp stars are 
modulated according to the visible magnetic field structure, point- 
ing to a defining role played by magnetic fields in exciting the os- 
cillations and shaping the main pulsation properties. Observation of 
the coincidence of the magnetic field and pulsation amplitude ex - 
trema gave rise to the oblique pulsator model (OPM. [Kurtzl l 1 982h . 
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According to this phenomenological framework, the main charac- 
teristics of non-radial pulsations in roAp stars can be attributed to 
an oblique axisymmetric dipole (£ = 1, m = 0) mode, aligned with 
the axis of a nearly axisymmetric quasi-dip olar magnetic field- 
Several theoretical inves t igations (e. g., Bigot & Dziembowski 
l2002t ISaio & Gautschvll2004 [Saioll2005h studied the effects of the 
distortion of oblique pulsation mode geometry by the global mag- 
netic field and stellar rotation. iBigot & Dziembowski ( 2002]) pre- 
dicted that the rotational distortion of pulsation eigenmodes is rep- 
resented by a superposition of spherical harmonic functions, con- 
taining substantial non-axisymmetric components, and that there is 
no alignment of th e pulsation axis and the d ipol ar magnetic field. 
On the other hand, ISaio & Gautschvl fc004h and lSaiol J2005T) dis- 
regarded effects of rotation, focusing on the detailed treatment of 
the interaction between p-modes and a strong magnetic field. These 
authors find an axisymmetric pulsation structure and alignment be- 
tween the magnetic field and the pulsation modes, and predict a 
strong confinement of pulsation amplitude to the magnetic field 
axis. 

Recently major progress in the observational study of roAp 
stars was achieved by employing high time resolution spec- 
troscopy in addition to the classical high-speed photometric mea- 
surements. Time-series spectroscopy of magnetic pulsators re- 
vealed a surprisi ng diversity in the pulsa t ional behaviour of dif- 
feren t lines (e.g., iKanaan &Hatzedll99l iKurtz. Elkin"&M athvs 
l200d and referen ces therein). Detailed analyses of the bright 
roAp star 7 Equ JSa vanov. Malanushenko & Ryabchikova 1 19991 ; 
Koch ukhov & Rv abchikova 2001a) demonstrated that spectro- 
scopic pulsational variability is dominated by the lines of rare-earth 
ions, especially those of Pr and Nd, which are strong and numerous 
in the roAp spectra. On the other hand, light and iron-peak ele- 
ments do not pulsate with amplitudes above 50-100 ms" 1 , which 
is at least an order of magnitude lower than the 1-5 km s~ 1 variabil- 
ity observed in the lines of rare-earth elements (REE). Many other 
roAp stars have been found to show a very similar overa l l pulsa- 
tional behaviour (e.g., Kochukhov & Rvabchikova! 1200 la iBalona 
l2002l;lMkrtichian. Hatzes & KanaarJI 20031) . 

The peculiar chara cteristics of the p-mode pu lsations in roAp 
stars were clarified by Rvabchi kova et al] d2002l) . who were the 
first to relate pulsational variability to vertical stratification of 
chemical elements. This study of the atmospheric properties of 
7 Equ showed that the light and iron-peak elements are en- 
hanced in the lower atmospheric layers, whereas REE ions are 
conce ntrated in a cloud with a lower boundary a t logrsooo < 
—4 |Mashonkina, Rvabchi kova & Rvabtsevl 120051) . Thus, high- 
amplitude pulsations observed in REE lines occur in the upper at- 
mosphere while lines of elements showing no significant variabil- 
ity form in the lower atmosphere. This leads to the following gen- 
eral picture of roAp pulsations: we observe a signature of a mag- 
netoacoustic wave, propagating outwards with increasing ampli- 
tude through the chemically stratified atmosphere. The presence of 
significant phase shifts bet ween the pulsation radial velocity ( RV) 
curves of different REEs JKochukhov & Rv abchikova 200 1_a|), o r 
even between lines of the same element l lMkrtichian et alj|2003l) . 
can be attributed to the chemical stratification effects and, possi- 
bly, to the short vertical wavelength of the running magnetoacoustic 
wave. These unique properties of roAp pulsations, combined with 
a presence of large vertical abundance gradients in the line-forming 
region, make it possible to resolve the vertical structure of p-modes 
and to study propagation of pulsation waves at the level of detail 
previously possible only for the Sun. 

The outstanding pulsational variability of REE lines in the 



spectra of rapidly rotating roAp stars permitted detailed mapping 
of the horizontal pulsation structure, leading to observational ev- 
idence that clearly favors one of the two aforementioned theoret- 
ical predictions. In fact, the short rotation periods of some roAp 
stars and the oblique nature of their non-radial oscillations al- 
low pulsational monitoring from different aspect angles, thus fa- 
cilitating reconstruction of the horizontal pulsation pattern. Using 
this unique geometrical property of roAp pulsations, and apply- 
ing pulsation Dopp l er im aging and moment analysis techniques, 
Kochukho\J j2004bU200d) has successfully disentangled different 
harmonic contributions to the pulsation geometry of the prototype 
roAp star HD 83368 and for the first time obtained an indepen- 
dent confirmation of the alignment of non-radial pulsations and 
magnetic field. It was shown that the mag netoacoustic pu l sation s 
in HD 83368 are shaped as suggested bv lSaio & Gautschvl d2004l) . 
whereas the non-axisymmetric pulsation c omponents predicted by 
the theory of Big ot & Dziembowskil ( 120021) cannot be detected. 

As for the slowly rotating roAp stars, despite dramatic re- 
cent progress in understanding the vertical structure of their pul- 
sation modes, relatively little attention has been paid to the prob- 
lem of inferring the horizontal geometry of pulsations. Typically 
it is assumed that a horizontal cross-section of non-radial pulsa- 
tion is given by an oblique axisymmetric mode of low angular de- 
gree, similar to the pulsation geometries inferred for rapidly ro- 
tating roAp stars, where periodic modulation of the geometrical 
aspect supplies additional information that helps to constrain the 
pulsation geometry. Thus, the question of systematic mode identi- 
fication, central to the studies of other types of pulsating stars, has 
not been thoroughly investigated in the case of sharp-lined mag- 
netic pulsators, which represent the majority of roAp stars. 

Understanding rapid line profile variations (LPV) of the 
slowly rotating roAp stars turns out to be a challenging task. 
The first observation that achi eved a signal-to-noise ratio sufficien t 
to detect profile variability ( Kochukhov & Rvabchi koval 1200 lal) 
demonstrated the presence of an unusual blue-to-red running fea- 
ture in the residual spectra of 7 Equ. Moreover, observations of a 
single-wave variability of the REE line width in 7 Equ is clearly in- 
consistent with any axisy mmetric pulsation geometry describ ed by 
sp herical harmonics (see | Aerts 1 De Pauw & W aelkens 1992). This 
led lKochukhov & Rvabchikovj ( 12001 al) to attribute LPV in 7 Equ 
to a non-axisym metric I = 2- 3 non-radial mode, contrary to the 
classical OPM of lKurtzUl982l) and in contrast to their own results 
for rapid rotators. Later Shibaha shi et alj d2004l) drew attention to 
the blue-to-red running wave in our time-resolved spectra of 7 Equ. 
These authors showed that the observed LPV is inconsistent with 
spectral variability expected for any low-degree mode in a slowly 
rotating star and, in turn, suggested an exotic shock wave model to 
explain the observed rapid line profile changes. 

In this study we embark on the task of detecting, characteriz- 
ing and interpreting LPV in sharp-lined roAp stars. Such a study 
has become feasible thanks to the availability of a large amount of 
high-quality, time-resolved spectroscopic data for a sizable sample 
of roAp stars. The aim of our analysis is to derive from observations 
a complete picture of LPV in a number of sharp-lined roAp stars 
and to tackle the problem of interpretation of the REE line variabil- 
ity. As a result of our roAp line profile variability survey, we reveal 
ubiquitous asymmetric LPV in many magnetic stars. We suggest 
a revised OPM that recognizes the influence of stellar pulsations 
on turbulence in the upper atmospheric layers of roAp stars, and in 
this way resolves the apparent discrepancy between the axisymmet- 
ric pulsation picture required by the standard OPM and asymmetric 
LPV observed in sharp-lined roAp stars. 
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Table 1. Fundamental parameters of program stars. 



HD 


Other 


T e ff 




v e sin i 


(B) 


P 


Reference 


number 


name 


(K) 




(kms" 1 ) 


(kG) 


(min) 




9289 


BW Cet 


7840 


4.15 


10.5 


2.0 


10.52 


Rvabchikova et al. (2007b) 


12932 


BN Cet 


7620 


4.15 


3.5 


1.7 


11.65 


Rvabchikova et al. (2007b) 


19918 


BTHyi 


8110 


4.34 


3.0 


1.6 


11.04 


Rvabchikova et al. (2007b) 


24712 


DOEri, HR 1217 


7250 


4.30 


5.6 


3.1 


6.13 


Rvabchikova et al. (1997) 


101065 


Przybylski's star 


6600 


4.20 


2.0 


2.3 


12.15 


Cowley et al. (2000) 


122970 


PPVir 


6930 


4.10 


4.5 


2.3 


11.19 


Rvabchikova et al. (2000) 


128898 


aCir 


7900 


4.20 


12.5 


1.5 


6.80 


Kupkaetal. (1996) 


134214 


HI Lib 


7315 


4.45 


2.0 


3.1 


5.69 


Ryabchikova et al. (2007b) 


137949 


33 Lib 


7550 


4.30 


^2.0 


5.0 


8.27 


Rvabchikova et al. (2004) 


201601 


7 Equ 


7700 


4.20 


<1.0 


4.1 


12.21-12.45 


Rvabchikova et al. (2002) 



The rest of the paper is organized as follows. Sect. [2] ex- 
plains our choice of the target roAp stars. Sect.[5]describes acqui- 
sition and reduction of time-series spectra. Investigations of spec- 
troscopic variability with residual spectra and moment analysis are 
presented in Sect. [4] Interpretation of the observations using theo- 
retical spectrum synthesis is given in Sect. [5] Results of our study 
are summarized and discussed in Sect. [6] 



2 TARGET SELECTION 

The pulsational LPV pattern of a roAp star strongly depends on the 
orientation of the pulsation axis with respect to the line of sight, 
and on the magnitude of Doppler shifts due to the stellar rota- 
tion. In the slowly rotating oblique pulsators (V p > v c sin i), the 
velocity field observed in REE lines is dominated by the pulsa- 
tional perturbations. These roAp stars are qualitatively similar to 
normal non-radially pulsating stars, except that the pulsation axis 
obliquity j3 plays the role of the inclination angle i. When the ro- 
tational Doppler shifts exceed the local pulsational velocity over a 
significant fraction of the visible stellar surface (V p <C v a sin i), the 
pulsational profile variability pattern changes its form, because it is 
now determined by a complex su perposition of the rotational veloc- 
ity field and oblique pulsations dKochuk hov 2004a). The behaviour 
of the line profil e moments is also modified in comparison to the 
sharp-lined stars (Kochukhov 2005). These circumstances suggest 
that different analysis strategies are best suited for fast and slow 
rotators. The former should be observed over the whole rotation 
period and modelled taking into account the complete oblique rota- 
tor/pulsator geometry. On the other hand, one can, at least as a first 
approximation, neglect rotation for the latter stars and treat them as 
normal pulsators. The scope of the present paper is observation and 
interpretation of the pulsational LPV in the second, slowly rotating, 
group of roAp stars. 

Based on our numerical simulations o f the line profile and 
moment variations jKochukhovll2004al l2005h . we have established 
v e sin i « 10-15 km s _1 (the exact figure depends on the pulsation 
geometry and amplitude) as an upper limit of the rotational Doppler 
shifts which can no longer be neglected in the spectroscopic mod- 
elling of oblique non-radial pulsations. This criterion, together with 
the availability of high-quality spectroscopic time-series data, has 
defined the selection of the program stars. Table Q] summarizes 
the main properties of ten roAp stars investigated in our paper. 
When available, we give T c g, logp, w e sini and the mean field 
modulus (B) derived in recent model atmosphere and chemical 
abundance analyses (see references in TableQJ. For 4 poorly stud- 



ied roAp stars (HD9289, HD 12932, HD 19918, and HD 134214) 
no detailed studies have been publ ished, hence we adopted at mo- 
spheric parameters determined by Rvabch ikova et al. 1 2007bl) us- 
ing Stromgren photometric calibrations jMoon & Dworetskvl 19851 : 
Napi wotzki. Schoenbemer & Wenskel 1 19931) . Estimates of rota- 
tional velocity and magnetic field for these four stars were obtained 
by fitting theoretical spectra to unblended spectral lines of differ- 
ent magnetic sensi t ivity. D etails of this procedure can be found in 
IRvabchikova et all d2007bh . 

Time-resolved spectra are also available for several roAp stars 
characterized by moderate to large line broadening: HD 42659, 
HD 60435, HD 80316, HD 83368, HD 84041, HD 99563. These 
stars do not satisfy our v c sin i criterion or have observations of 
low quality (HD 60435), and were not included in the analysis. 

We note that some objects selected according to their low 
v e sin i may be relatively rapid rotators seen close to one of the 
oles. The 3.877 d rotation period obtained by Rvab chikova et al] 
20051) for HD 12 2970 and the 4.471 d rotational modulation found 
for H P 128898 dKurtz et alj 1 19941 : iKochukhov & Rvabchikov3 
2001b) hint on such a viewing geometry. Thus, an additional un- 
certainty may enter interpretation of the spectroscopic variability 
of HD 122970, HD 128898 and HD9289 (whose rotation period 
is unknown, but cannot be long, given a v e sin i of 10.5 kms -1 .) 
Nevertheless, spectroscopic observables studied in our paper are 
mostly sensitive to the projected rotational velocity and should 
not be strongly distorted by the special orientation of HD 122970, 
HD 128898 and similar objects. Consequently, we make no explicit 
distinction between intrinsically slow rotators and pole-on sharp- 
lined stars. 

The pulsation periods listed in Table[T]correspond to variabil- 
ity due to the dominant pulsation mode at the time of observations. 
In this study we do not perform a detailed frequency analysis, and 
only deal with short (typically ~ 2h-long) spectroscopic time se- 
ries. Therefore, possible small inaccuracies in the adopted pulsation 
periods have no influence on the results presented below. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Echelle spectra 

The main spectroscopic data set analysed in our study includes 958 
observations of 8 roAp stars obtained with the UVES instrument 
at the ESO VLT between October 8, 2003 and Mar ch 12, 2004 
in the context of the observing program 072.D-0138 dKurtz et al] 
2006). The ESO Archive facility was used to search and retrieve 
science exposures and the respective calibration frames. For these 
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Table 2. Journal of time-series spectroscopic observations of roAp stars. Superscript numbers 1-6 mark different data sets of 
HD 201601. 



HD 


Instrument/ 


Stall HJD 


End HJD 


Wavelength 


Number of 


Exposure 


Overhead 


Peak 


number 


telescope 


(2450000+) 


(2450000+) 


region (A) 


exposures 


time (s) 


time (s) 


SNR 


9289 


UVES/VLT 


2920.54506 


2920.62881 


4960-6990 


111 


40 


25 


90 


12932 


UVES/VLT 


2921.62234 


2921.70532 


4960-6990 


69 


80 


25 


90 


19918 


UVES/VLT 


2921.52607 


2921.60905 


4960-6990 


69 


80 


25 


100 


24712 


UVES/VLT 


3321.65732 


3321.74421 


4960-6990 


92 


50 


22 


300 


101065 


UVES/VLT 


3071.67758 


3071.76032 


4960-6990 


111 


40 


25 


180 


122970 


UVES/VLT 


3069.70977 


3069.79359 


4960-6990 


111 


40 


25 


160 


128898 


UVES/VLT 


3073.80059 


3073.88262 


4960-6990 


265 


1.5 


25 


250 


134214 


UVES/VLT 


3070.77571 


3070.85848 


4960-6990 


111 


40 


25 


260 


137949 


UVES/VLT 


3071.76312 


3071.84598 


4960-6990 


111 


40 


25 


350 


201601 


CES/ESO 3.6-m 1 


1381.78344 


1381.82390 


6140-6165 


31 


60 


52 


190 


201601 


Gecko/CFHT 2 


2186.70618 


2186.80296 


6543-6656 


64 


90 


43 


250 


201601 


Gecko/CFHT 2 


2186.82456 


2186.92308 


6619-6729 


65 


90 


43 


230 


201601 


Gecko/CFHT 3 


2540.72237 


2540.83199 


5283-5343 


65 


90 


53 


160 


201601 


Gecko/CFHT 3 


2540.86101 


2540.96843 


5821-5887 


65 


90 


53 


130 


201601 


Gecko/CFHT 4 


2541.85433 


2541.93164 


5821-5887 


50 


90 


43 


140 


201601 


Gecko/CFHT 5 


2542.82688 


2542.95051 


6542-6657 


60 


90 


43 


170 


201601 


Gecko/CFHT 6 


2543.82314 


2543.92191 


6104-6194 


64 


90 


43 


230 



observations the red arm of the UVES spectrometer was configured 
to observe the spectral region 4960-6990 A (central wavelength 
6000 A). The wavelength coverage is complete, except for a 100 A 
gap centred at 6000 A. Observations were obtained with the high- 
resolution UVES image slicer (slicer No. 3), providing improved 
radial velocity stability and giving a maximum resolving power of 
A/AA ~ 115 000. 

Observations of each target cover 2 hours and consist of an 
uninterrupted spectroscopic time series with a total number of ex- 
posures ranging from 69 to 265. The length of individual exposures 
was 40 s or 80 s, except for the brightest roAp star HD 128898 
(a Cir), for which a 1.5 s exposure time was used. The ultra-fast 
(625kHz/4pt) readout mode of the UVES CCDs made it possible 
to limit overhead to ~ 25 s, thus giving a duty cycle of 70-80% 
for the majority of targets. The signal-to-noise ratio of individual 
spectra is between 90 and 350, as estimated from the dispersion 
of the stellar fluxes in the line-free regions. Detailed description of 
observations for each star is given in Table[2] 

The red 600 nm UVES data set was complemented by the 
observations of HD 24712 obtained on November 12, 2004 on 
DDT program 274.D-501 1. 92 time-resolved spectra were acquired 
with the UVES spectrometer, which was used in the 390+580 nm 
dichroic mode (wavelength coverage 3300-4420 A and 4790- 
6750 A). A detailed de scription of the acquisitio n and reduction 
of these data is given by Rvabchikova et al. ( 2007a). 

All echelle spectra were reduced and normalized to the contin- 
uum level with a routine specially developed by one of us (DL) for 
fast reduction of spectroscopic time-series observations. A special 
modificatio n of the Vienna automatic p ipeline for echelle spectra 
processing iTsvmbalT Lvashko & Weiss 2003) was developed. All 
bias and flat field images were median averaged before calibration. 
The scattered light was subtracted using a 2-D background approxi- 
mation. For cleaning cosmic ray hits we applied an algorithm which 
compares the direct and reversed spectral profiles. To determine the 
boundaries of echelle orders, the code used a special template for 
each order position in each row across the dispersion axis. The shift 
of the row spectra relative to the template was derived by a cross- 
correlation technique. Wavelength calibration for each roAp star 
was based on a single ThAr exposure, recorded immediately after 



the respective stellar time series. Calibration was done by a 2-D 
approximation of the dispersion surface. An internal accuracy of 
30-40 m s _1 was achieved by using several hundred ThAr lines in 
all echelle orders. The final step of continuum normalization and 
merging of echelle orders was carried out by transformation of the 
flat field blaze function to the response function in each order. 

The global continuum normalization was improved by itera- 
tive fitting of a smoothing spline function to the high points in the 
average spectrum of each star. With this procedure we have cor- 
rected an underestimate of the continuum level, unavoidable in the 
analysis of small spectral regions of the crowded spectra of cool 
Ap stars. Correct determination of the absolute continuum level 
is important for retrieving unbiased amplitudes of radial velocity 
variability and for studying LPV. In addition to the global contin- 
uum correction, the spectroscopic time series were post-processed 
to ensure homogeneity in the continuum normalization of individ- 
ual spectra. Each extracted spectrum was divided by the mean spec- 
trum. The resultant ratio was heavily smoothed, and was then used 
to correct the continuum in the individual spectra. 

3.2 Single order spectra 

The time-resolved observations of HD 201601 were obtained in 
2001-2002 using the single-order //4 Gecko coude spectrograph 
with the EEV1 CCD at the 3.6-m Canada-France-Hawaii telescope. 
In total we gathered 491 observations, covering several spectral 
windows (see Tabled in the 5823-6729 A wavelength region. Ob- 
servations have a resolving power, determined from the widths of a 
number of ThAr comparison lines, of about 115 000. In all obser- 
vations of HD 201601 we used an exposure time of 90 s. Individual 
exposures have signal-to-noise ratio between 130 and 230. 

The spectra were reduced using standard IRAF tasks. Each 
stellar, flat and calibration frame had a mean bias subtracted and 
was then cleaned of cosmic ray hits and extracted to one dimen- 
sion. Extracted stellar spectra were divided by an extracted mean 
flat field, and the continuum was fit with a third-order Legendre 
polynomial, using the same rejection parameters for all spectra so 
that the continuum fit is as uniform as possible. The wavelength 
scale was established using about 40 lines of a ThAr emission lamp, 
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resulting in an RMS scatter about the adopted pixel-wavelength 
polynomial (a sixth-order Legendre polynomial) of about 30 m s~ 1 . 
The wavelength scale was linearly interpolated between ThAr lamp 
spectra taken before and after the stellar series. The spectra were 
not resampled to a linear wavelength spacing. 

In this study we also used 31 very high resolution time- 
resolved observations of HP 20 1601 (7 Equ) analysed by 
iKochukhov & Rvabchikovi feOOld) . These data were obtained in 
1999 with the Coude Echelle Spectrograph (CES), fiber-linked to 
the Cassegrain focus of the ESO 3.6-m telescope. The highest reso- 
lution CES image slicer and the ESO CCD#38 were used, allowing 
us to reach a resolving power of A/ A A = 166 000 and record spec- 
tra in the 6140-6165 A wavelengt h interval. We refer the reader to 
IKochu khov & Rvabc hikoval j2001ah for other details of the acqui- 
sition and reduction of the CES spectra of HD 20 1601. 

Post-processing of the extracted spectra of HD 201601 was 
done consistently with the procedure adopted for UVES observa- 
tions. 



4 INVESTIGATION OF SPECTROSCOPIC 
VARIABILITY 

4.1 Choice of spectral lines 

Non-radial pulsations in roAp stars are best observed using lines 
of the rare-earth ions. Here we aim at detailed investigation of the 
pulsational LPV and, therefore, are limited to the analysis of vari- 
ability in a few of the strongest REE lines which show the clear- 
est signatures of oscillations. Our experience, and the results of 
previous studies, suggest that strong lines of Pr III and Nd III are 
most useful for the detection and interpretation of the rapid LPV. 
Therefore, we have chosen a sample of unblended doubly ionized 
lines of Pr and Nd. In addition, we studied a few weaker lines 
of Nd II and Tb III in order to sample different layers in the up- 
per atmospheres of roAp stars a nd to resolve pos sible depth de- 
pendence of LPV. The VALD feupka et alj|l999h and DREAM 
(Biemont, Palmeri & Ouine databases, complem ented with 

the extended Nd III line list ( Rvabchikova et al. 2006), were used 
as a basis for line identification. In total, we looked at 17 REE 
lines, including Prill 5284, 5300, 6160, 6706 A, Ndn 5311, 5319, 
6650 A, Ndm 5286, 5294, 5802, 5851, 6145, 6550, 6691 A, and 
Tbm 5505, 5847, 6323 A. Not all lines can be studied in each 
star because of blending considerations, different spectral cover- 
age of the time-resolved spectra, and the intrinsic star-to-star vari- 
ation in pulsational behaviour. We preferred to analyse lines with 
A ^ 6000 A in stars observed with UVES, whereas red REE lines 
with central wavelengths above 6000 A were studied in the CFHT 
observations of 7 Equ. In fact, an exhaustive investigation of all 
the lines listed above in each star is unnecessary because the infor- 
mation content of many REE features is redundant (this especially 
concerns the numerous strong Nd III lines). 

A remarkable diversity of the pulsational amplitudes and 
phases observed in the spectral lines of different chemical ele- 
ments indicate the importance of the depth dependence of chem- 
ical abundances and of pulsation wave properties in roAp atmo- 
spheres. Taking this into account, we expect to find a significant 
depth dependence of the pulsational LPV. This is why it is im- 
portant to study variability in diverse REE lines, and to associate 
observed LPV patterns with specific atmospheric layers. However, 
the question of the formation depth of REE lines in the atmo- 
spheres of roAp stars is an extremely complex one. Ideally it should 



be ad dressed with detailed NL TE chemical stratification calcula- 
tions (Mas honkina et al]|2005h . possibly taking into account the 
anomalous structure of the upper atmo spheres of cool Ap stars 
( IKochukhov. Bagnulo & Barkleml 120021) . Such an analysis is ex- 
tremely demanding in terms of computing time and input physics, 
and up to now has been carr ied out only for Nd II a nd Nd III in two 
representative roAp stars jMashonkina et alj|2005l) . The results of 
these elaborate NLTE-chemical stratification calculations show that 
line intensity can be used as a good proxy for the relative forma- 
tion depth of spectral features belonging to the same element. For 
instance, weaker singly ionized Nd lines are formed deeper than 
strong Nd III lines. On the other hand, the intensity of the lines of 
different elements is not a reliable indicator of the formation depth 
because there is no reason to assume that vertical stratification is 
the same for all REE species. Nevertheless, one can still assign 
an approximate relative formation depth to the lines of different 
REEs based on phases of RV variability. In the framework of a run- 
ning magnetoacoustic wave which propagates outwards, the lines 
showing later RV maximum are formed higher in the atmosphere. 
We will use this property for intercomparison of the results ob- 
tained for Pr, Nd and Tb lines. One has to keep in mind, however, 
that different horizontal sampling of the pulsation velocity field 
due to dissimilar spotted distribution of different REEs may distort 
this phase-depth relation. Moreover, any judgment about formation 
heights based solely on the time of the velocity extrema becomes 
meaningless when the vertical p-mode structure is dominated by 
the standing wave component, and pulsation nodes with the asso- 
ciated 180° phase change are observed in the line-forming region. 
One roAp star includ ed in our sample, HD 137949, is known to 
show this behaviour (Mkrtic hian et alj2 003). 

4.2 Residual spectra 

We start analysis of the pulsational LPV in roAp stars by investigat- 
ing the behaviour of the residual spectra as a function of pulsation 
phase. This is the most general approach to characterize LPV of 
pulsating stars and, because the whole information content of the 
line profile changes is considered, this method may be preferable 
in some cases to analysis of average secondary quantities, such as 
line profile moments and bisectors. On the other hand, analysis of 
the residual spectra has a drawback, in that it is impossible to di- 
rectly relate features in the LPV pattern to physical properties of 
non-radial pulsations. LPV represents a result of a non-trivial con- 
volution of the intrinsic stellar line profiles and velocities due to 
pulsation and stellar rotation. Therefore, credible interpretation of 
the residual profile variation is possible only with the help of spec- 
trum synthesis calculations which adequately describe the afore- 
mentioned effectfl 

Investigation of the residual spectra is feasible only for time- 
resolved observations of sufficient quality. The signal-to-noise ra- 
tio of the time-series spectra of roAp stars available to us is quite 
uneven, which precludes direct detection of the periodic LPV in 
fainter stars (about half of our sample). To circumvent this diffi- 
culty, we developed a special spectral coaddition procedure. First, 
we derive an average spectrum for each star and subtract it from the 
individual observations. Then pulsation phases are computed using 

1 In this context we note that the assertion by Shibahashi et al. 
120041) that certain featu r es of the LPV observed in 7 Equ by 
IKochukhov & Rvab chikova 1 2001a) can be interpreted in terms of the pul- 
sation velocity amplitude is erroneous. 
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the periods listed in Table [T] Observations with similar phases are 
averaged. We typically divided the pulsation cycle into 10 phase 
bins, except for the CFHT and ESO 3.6-m 7 Equ data sets where 
1 1 and 13 phase bins were used, respectively. 

Variation of the residual profiles of the Ndm 6145 A and 
Pr III 6160 A lines in 7 Equ, processed with our phase binning pro- 
cedure, is illustrated in Fig. Q] In this and other similar plots we 
show the average line profile on top, the time series of the differ- 
ence spectra is plotted in the middle (with time increasing down- 
ward), and the standard deviation profile is shown at the bottom 
of the plot. Comparison of the binned LPV presented in Fig. [Tj 
w ith the profile variations in the ori ginal observations (see fig. 1 
in lKochukhov & R vabchikova 2001a) shows very good agreement. 
This indicates that application of the phase binning procedure to 
the time-resolved spectra of roAp stars is able to enhance S/N 
while retaining the original LPV pattern. As was observed by 
iKochukhov & Rvabchikovj feOOlal) . we find that LPV of the Nd III 
and Pr III lines studied is asymmetric, and is dominated by the blue- 
to-red moving feature. 

The LPV pattern detected in the CFHT observations of 
7 Equ is shown in Fig. [2] This 7 Equ data set is unique 
in that it shows an exceptionally high amplitude of the pulsa- 
tional variability (RV semi-amplitude of up to 1 kms -1 for a 
number of REE lines) compar ed to all other observations of 
7 Equ reported in the literature (K ochukhov & Rvabchikova 2001a; 
IKochukhov. Rvabchikova & Piskunovll2004 ISavanov et alj|2005h . 
Fig.[2]compares behaviour of the core of Ha, Prill 6706 A, Ndm 
6550, 6691 A and Ndll 6650 A. (Two plots are shown for the lat- 
ter line, corresponding to observations in the two consecutive time 
series recorded during the same night.) One can see a clear differ- 
ence in the LPV pattern displayed by the Ha core, Nd II 6650 A, 
and Nd III 6691 A on the one hand and by the Nd III 6550 A, Prill 
6706 A features on the other hand. The first group of lines shows 
nearly symmetric, S-shaped variability, whereas the two lines from 
the second group behave similarly to the Ndm 6145 A and Prill 
6160 A lines in Fig. [T] and show a blue-to-red moving wave. 
We suggest that this difference reflects the depth variation of the 
pulsational LPV pattern. Sophisticated line formation calculations 
(Mashonkina et al. 1 120051) predict that the absorption in the Ndll 
6650 A line originates at the bottom of the Nd-rich layer. The Nd III 
669 1 A line is one of the few doubly ionized Nd lines whose forma- 
tion depth overlaps with that of Nd II. This is confirmed by a rather 
small delay (Aip = 0.08, hereafter we give phases and phase dif- 
ferences in units of pulsation period) of its RV maximum relative to 
the radial velocity variation of Nd II 6650 A. In contrast, the strong 
line Ndm 6550 A forms substantially higher and shows a large 
phase delay (Aip = 0.21). Detailed NLTE radiative transfer calcu- 
lations are not yet available for Pr III, but, given that the weak Pr III 
6706 A line has a phase lag similar to that of Nd III, its effective for- 
mation height should not differ much. Indeed, we see qualitatively 
the same LPV pattern for this Pr III feature and for Nd III 6550 A. 

The formation depth of the Ha core is difficult to estimate. 
Based on its symmetric LPV pattern, we expect the Ha core to 
form below the Nd-rich layer, in agreement with the proposed depth 
dependence of LPV in metal lines. Nevertheless, the large intrinsic 
Doppler width of the Ha core does not permit a straightforward 
comparison with the variability detected in REE lines. In addition, 
the pulsational changes of the Ha core are too low in some of the 
roAp stars included in our study. For these reasons we refrain from 
detailed analysis of the Ha profile behaviour in this paper. 

The CFHT observations of 7 Equ discussed in this section 
constitute the first detection of the depth variation of the pulsa- 
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Figure 1. Profile variations of the Ndm 6145 A and Pr III 6160 A lines 
in the CES observations of 7 Equ obtained on 22 June, 1999. The average 
spectrum is plotted in the upper part of each panel. The difference spectra 
as a function of pulsation phase are presented in the middle part. Observa- 
tions for consecutive pulsation phases are shifted in the vertical direction, 
and each phase is shown twice. The bottom curve presents the wavelength 
dependence of the the standard deviation (expanded by a factor of 10). 



tional LPV pattern in a roAp star. Moreover, to our knowledge, 
such behaviour has not previously been observed in any main se- 
quence non-radially pulsating star. We conclude that the pulsational 
LPV changes significantly within the REE-enriched cloud in the 
upper atmosphere of 7 Equ: the symmetric, S-shaped LPV pattern 
of the lines from deeper layers is replaced by asymmetric, blue-to- 
red moving features in the LPV of lines formed higher up. 

For the analysis of the REE line profile variability in our main 
UVES data set of nine roAp stars we used the Nd II 5319 A, Nd III 
5294 A, Pr III 5300 A and Tb III 5847 A lines. The Nd II feature is 
sufficiently strong to be observed in all stars and is formed at the 
bottom of the Nd-rich cloud. The 5294 A line is one of the strongest 
Nd III lines in the optical spectra of cool Ap stars. This line samples 
high atmospheric layers. In general, nothing can be said about the 
formation heights of the chosen lines of doubly ionized Pr and Tb 
because they depend on currently unknown details of stratification 
of these REE species in individual stars. 

Fig. [3] shows residual profile variation for the Ndll 5319 A 
line. Among the nine roAp stars presented in this plot, seven show 
a symmetric or nearly symmetric (HD 19918), S-shaped LPV pat- 
tern similar to the one which we found for Ndll lines in 7 Equ. 
HD9289 shows a weak red-to-blue moving wave - an odd be- 
haviour which is not seen in any other star. Time series of the resid- 
ual spectra are inconclusive for HD 137949 because at the heights 
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Figure 2. Profile variations of the Ndm 6550, Ha core, Ndll 6650, Ndm 6691, and Prill 6707 A lines in the CFHT observations of 7 Equ obtained on 4 
October, 2001. The layout of the figure is similar to Fig.[T] 
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Figure 3. Profile variations of the Ndll 5319 A line for nine roAp stars. The layout of the figure is similar to Fig.[T] except that here an arbitrary intensity scale 
is used for the residual and standard deviation spectra. The vertical dashed lines indicate position of line centres. 



sampled by Ndll 5319 A, a pulsation node is observed in this star 
and the variation is dominated by the first harmonic of the main fre- 
quency. How^vet^_die_Ha core in HD 137949 has symmetric vari- 
ability dKurtz et al.l2004l) . quite similar to the behaviour of Nd II in 
the majority of the roAp stars studied. 

Considering variation of the Ndm 5294 A line (Fig. 
we detect a well-developed LPV pattern with blue-to-red moving 
waves in HD9289, HD 12932, HD 19918, HD 24712, HD 128898, 
HD 134214 and HD 137949. The phase lag of this Nd III line with 
respect to the RV changes of Nd II 5319 A is typically in the range 
0.17-0.22, but smaller (0.15-0.17) for HD24712 and HD 134214. 
These two stars exhibit LPV of a transitional type: it is possi- 
ble to trace a blue-to-red moving wave; however, the latter is 



not continuous, but formed by two separate 'bumps' in the time- 
sequence of the residual spectra. The two cooler stars, HD 101065 
and HD 122970, retain the S-shaped LPV which was present in 
the Ndll 5319 A line. Remarkably, these two stars also show a 
much smaller phase lag (Atp = 0.01-0.04) between Nd III and Nd II 
compared to the rest of the roAp star sample. Either the forma- 
tion heights of the Nd II and Nd III lines are less different in cooler 
roAp stars due to a special form of Nd stratification, or the pulsation 
wave structure sampled by the Nd lines is dominated by a standing 
wave. It is worth m entioning that abundance analysis of HD 101065 
dCowlev et alj2000l) provides no evidence of the Nd II-Nd III abun- 
dance anomaly found in other roAp stars. 

The variation of the Prill 5300 A line is illustrated in Fig. [5] 
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Figure 4. Profile variations of the Nd III 5294 A line for nine roAp stars. The layout of the figure is similar to Fig. [5] 
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Figure 5. Profile variations of the Pr III 5300 A line for nine roAp stars. The layout of the figure is similar to Fig. [5] 



The negligible phase lag of this REE line with respect to Ndm 
5294 A in HD9289, HD 12932 and HD 19918 corresponds to es- 
sentially the same LPV behaviour as shown by the Ndlll line. 
Thus, we may suggest that Prill and Nd III lines in these stars form 
in similar layers - a situation reminiscent of the PrIIl/Ndlll line 
variation in 7 Equ. HD 24712, HD 128898 and HD 134214 show 
a larger phase lag (Aip = 0.06-0. 11) and a stronger asymmetry in 
the Prill line variability. HD 101065 and HD 122970 continue to 
show symmetric LPV pattern. Analysis of the Prill 5300 A line 
in HD 137949 is complicated by the low RV amplitude and by the 
presence of harmonic oscillations. Nevertheless, blue-to-red mov- 
ing features are still visible in the time evolution of the residual 
spectra constructed for this star. 

Finally, we investigated variation of the Tbm 5847 A line. 
Its pulsational behaviour is diverse in the studied sample of roAp 
stars. Large pulsational amplitude and significant phase shifts of 
this line in some stars may indicate that Tb III probes atmospheric 
layers above those to which Ndlll and Prill lines are sensitive. 



Fig. [6] shows variation of Tbm 5847 A in four cooler roAp 
stars (HD24712, HD 101065, HD 122970, HD 134214). A very 
clear asymmetric LPV and significant (0.16-0.23) phase delay 
with respect to Prill 5300 A is found for HD 24712, HD 101065 
and HD 134214. Variation of the Tbm line in HD 122970 ex- 
hibits a transitional behaviour (cf. the Ndlll 5294 A variation in 
HD 134214), showing the first signature of the blue-to-red running 
waves. 

Our survey of the pulsational LPV in Nd, Pr and Tb lines in 
nine roAp stars reveals ubiquitous signatures of the depth variation 
of the profile variability pattern in every star. This corroborates re- 
sults obtained for 7 Equ and shows that a large vertical gradient of 
the amplitude and phase of pulsation waves observed in the atmo- 
spheres of roAp stars is associated with a similar rapid change of 
the LPV pattern. All stars, except HD 9289, show symmetric LPV 
in deeper layers. All stars, without exception, ultimately develop 
asymmetric LPV in the upper atmospheric region. This behaviour is 
apparently correlated with T e g : we find that hotter stars (HD 9289, 
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Figure 6. Profile variations of the Tb III 5847 A line for four cool roAp 
stars. The layout of the figure is similar to Fig. [5] 

HD 12932, HD 19918, HD 128898, HD 137949) develop interest- 
ing asymmetric LPV pattern lower (in the layers sampled by Nd III 
5294 A) than cooler stars (HD 24712, HD 101065, HD 122970, 
HD 134214). We also identified a characteristic change in the shape 
of the standard deviation profile associated with the transition from 
symmetric to asymmetric LPV. In the former case the blue and 
red peaks in the standard deviation spectra are well separated and 
there is almost no variation in the line core. In the latter case, 
the two peaks are often merged together and some residual inten- 
sity changes are observed in the line centre (e.g., Prill 5300 A 
in HD 19918 and HD 134214, TbiII 5847 A in HD24712 and 
HD 134214). 

4.3 Line profile moments 

Interpretation of the line profile moment variation of slowly rotat- 
ing pulsating s t ars is a useful puls a tion mode diagnostic te chnique 
dBalonal 1 19861 ; lAerts et alj 1 19921; iBriquet & Aertsl 120031) . Many 
studies have successfully utilized analysis of the pulsational be- 
haviour of the low-order line profile moments to determine os- 
cillation frequenc ies and to identify pulsation modes. Recently 
lKochukhovl (2005) has generalized the moment method to treat pul- 
sations arbitrarily inclined with respect to the stellar rotation axis. 
This modification of the moment technique is essential for its ap- 
plication to roAp stars. 

In the present study we complement direct analysis of the pul- 
sational LPV based on residual spectra with the study of pulsational 
perturbation of the REE line profile momen ts. These quantiti es are 
defined according to the following relation dAertsetal] 19921) 

{yJ) = ^ f^MY (i - s0 /£a - Si) , (i) 

where (V 1 ) is the line profile moment of jth order, A^ and Si 
are the wavelength and normalized intensity at the ith pixel of ob- 



served spectrum and c is the speed of light. Moments are computed 
with respect to the wavelength (A) corresponding to the centre- 
of-gravity of absorption feature in the average spectrum. The first 
moment (V) represents a measure of radial velocity, whereas the 
second moment (V 2 ) is related to line width. The sum in Eq. QJ is 
cut off when intensity reaches continuum level or at the point when 
blending by the neighboring lines becomes significant. To improve 
the accuracy of the calculation of moments with Eq. Q}, we have 
resampled the spectra of roAp stars into a fine wavelength grid with 
the help of cubic splines. In each star the pulsational variation of the 
first and second moments was investigated for a set of unblended 
lines of Pr III, Nd II, Nd III and Tb III. 

Time series of the line profile moments are analysed with the 
least-squares fitting method in order to characterize variation at the 
principal frequency and its first harmonic, and to extract the re- 
spective pulsation phases. Moment variation is approximated by 
the expression 

2 

ly 3 ) = Vo(i-io) +^2VjiCos{2Tr[i(t-t )/P-(pji]}. (2) 

i = l 

Here the first term takes into account possible drift of the spec- 
trograph's zero point. The HJD of the first exposure of the star at 
a given night (see Table [2]l was chosen as the reference time to- 
Vji and tfiji are, respectively, the amplitude and phase of the jth 
moment variation with the main period (i = 1) and half the main 
period (i = 2). With the minus sign in front of ipji a larger phase 
corresponds to a later time of maximum of (V) or (V 2 ). This phase 
convention is natural when discussing effects of the outward prop- 
agation of pulsation waves in the atmospheres of roAp stars. The 
pulsations periods in Eq. {2} are fixed to the va lues listed in Table[T] 

From theoretical considerations jAerts et al.l 1 19921 : 
iKochukhovl l2005h it is expected that the RV of the star ex- 
hibiting low-amplitude linear non-radial pulsations changes 
sinusoidally, whereas variation of the second moment may contain 
a significant contribution of the first harmonic. The amplitude ratio 
V21/V22 characterizing the overall shape of the (V 2 ) phase curve 
is very sensitive to the horizontal structure of the pulsation velocity 
field. However, possible non-linearities may lead to a non-zero 
harmonic contribution to the first moment variation and could 
distort the second moment amplitude ratio diagnostics. Thus, in 
order to assess possible contribution of the intrinsic harmonic 
variability, we fitted the same expression given by Eq. (O to the 
time series of both (V) and (V^ 2 ). 

An example of the moment analysis for several REE lines in 
HD 19918 and HD 101065 is displayed in Fig.[7] The first and sec- 
ond moment measurements are compared with the best-fitting co- 
sine curves computed according to Eq. {2j. Numerical results of 
fitting the first and second moment variability are compiled in Ta- 
bles[3]and|4]for the roAp stars observed with UVES and in Table|5] 
for HD 201601. We have verified that the pulsational behaviour of 
the latter star does not change significantly in the course of individ- 
ual observing nights. Thus, we consider moment measurements ob- 
tained in different wavelength regions observed on the same nights 
to be part of the same data set. This gives us six groups of time- 
resolved spectra, referred to by the numbers 1-6 in Tables [2] and 

s 

In the remaining part of this section we explore various trends 
in the amplitudes and phases of moment variation and relate these 
quantities to the LPV of REE lines discussed in Sect. 14.21 We 
start by noting that majority of the roAp stars in our sample 
show nearly sinusoidal RV variation. Marginally significant non- 
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Figure 7. Variation of line profile moments for rare-earth lines in a) HD 19918 and b) HD 101065. For each star the behaviour of the first (V) (left panels) 
and second (V 2 ) (right panels) profile moments is presented. Symbols show the moment measurements as a function of pulsation phase. The curves illustrate 
a fit including variation with the principal frequency and its first harmonic. 



zero V12 amplitudes are obtained only for a few REE lines in 
HD9289, HD 101065, HD 122970 and HD 128898. Nevertheless, 
for all these objects condition Vu 3> V12 is fulfilled, indicating 
that possible non-linearity or a contribution arising from the res- 
onantly excited harmonic oscillation is negligible. Thus, the pres- 
ence of harmonic variation of the second moment and the amplitude 
ratio V21 / V22 retains their value as mode structure diagnostics. An 
atypical behaviour is found for HD 137949. A sizable fraction of 
REE lines in this star shows significant, in some cases dominant, 
double-wave RV variation and correspondingly have large V12 am- 
plitude. In this situation, the origin of harmonic variability of {V 2 } 
is ambiguous. 

Following arguments presented in Sect. 14. II we use the phase 
<pn of RV fluctuations as a proxy of the relative formation height 
for REE lines. For pulsation waves propagating outwards in the 
radial direction, a larger RV phase corresponds to a later time of 
maximum velocity, which suggests a higher physical location of 
the line-forming region of the REE ion. 

For all stars except HD 137949 we discover a dramatic change 
of the second moment amplitude V21 with height. This trend is ev- 
ident from inspection of Fig.|7]for HD 19918 and HD 101065. The 
Nd II lines from deeper layers often exhibit large, or at least measur- 
able harmonic variation of {V 2 ). As the pulsation wave propagates 
towards higher layers, the V21 amplitude increases steeply and the 
shape of {V 2 ) phase curve transforms from a low-amplitude, pre- 
dominantly double wave to a large-amplitude, mostly single wave. 
This tendency is summarized for all stars in Figs. [8] and [9] We 
find a persistent increase of the principal frequency component 
in the (V 2 ) variation with increasing phase (pn. Typically, Ndll 
lines have low V21 amplitude; this parameter increases in strong 
Nd III lines, continues increasing for Pr III and finally reaches max- 
imum in Tbm (HD 12932, HD24712, HD 101065, HD 122970, 
HD 128898, HD 134214). For a few stars (HD9289, HD 19918 and 
HD 201601) where TbiII lines do not show large-amplitude vari- 
ability of {V 2 }, the V21 amplitude maximum occurs in Prill lines. 

Comparing the inferred depth dependence of the second mo- 
ment variation with the line-to-line changes of LPV pattern, we 
discovered a remarkable correspondence between high-amplitude, 
single- wave pulsational variability of (V 2 ) and the presence of a 



blue-to-red wave pattern in the LPV. Detailed intercomparison of 
the residual profile time series displayed in Figs. [3}j6] and the mo- 
ment fitting results reveals that all those lines, in particular Pr III 
and Tb III, showing asymmetric LPV are invariably characterized 
by a large V21 amplitude and high V21 jV\\ ratio. The latter param- 
eter reaches a value of 10-15 for Pr III and Tb III in most stars. The 
depth variation of LPV pattern and V21 amplitude is also consistent. 
The two stars, HD 101065 and HD 122970, that develop asymmet- 
ric LPV late and show it prominently only in the Tb III lines also ex- 
hibit less persistent single-wave second moment changes and have 
lower V2\jV\\ ratios (2-7). Thus, an important conclusion emerg- 
ing from our analysis of the first and second moment variation with 
pulsation phase is that a large-amplitude, single-wave behaviour of 
(V^ 2 ) and asymmetric LPV pattern are closely related and proba- 
bly have the same physical origin. We stress that this correlation 
between the two distinct characteristics of the roAp spectroscopic 
variability is universal and is found in all stars where intrinsic har- 
monic variability is low enough not to distort the moment tech- 
nique diagnosis. Consequently, any theory purporting to explain 
the asymmetric LPV pattern of sharp-lined roAp stars must also 
reproduce the main features of the moment behaviour. 

Further insight into pulsational line width perturbations is ob- 
tained by considering the phase shift between the dominant single- 
wave variability of (V) and { V 2 ) . In Figs.[8]and[9]we plot the phase 
delay, (^21 — ipu, of the second moment with respect to the phase 
tpii of the RV curve. The value of this phase lag is diverse for 
different groups of spectral lines, and no consistent depth depen- 
dence is found. However, an interesting coherency is found for the 
numerical value of tp2i — ipu in lines with dominant single-wave 
variation of the second moment. Using lines of Pr III and Tb III (for 
HD 201601 the Nd III lines are also considered) which are sensitive 
to pulsations in the outermost stellar layers, we find that ip2i — tpii 
is always positive (see average values of the phase lag listed in Ta- 
bles[2j-[5]>. This means that the (V 2 ) maximum occurs later in time 
compared to the RV maximum. Furthermore, the ip2i — fin param- 
eter spans a small range of values between 0.19 and 0.28. (The few 
exceptions include HD9289, which has (ip2i — fin) = 0.09, and 
one 7 Equ data set with (9521 — fin) = 0.15). For seven out of 
ten stars studied, (tp2i — 'fin) is consistent within la confidence 
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Figure 8. The amplitude V21 of the second moment variation at the main 
frequency (left panels) and the phase difference ip2i — Pu (right pan- 
els) as a function of the pulsation phase ipn of RV changes for HD 9289, 
HD 12932, HD 19918, HD 24712 and HD 101065. Symbols correspond to 
different REE ions: Prill (A), Ndll (❖), Ndlll (•), Tbm (*). 



level with exactly one-quarter period phase lag between RV and 
(V 2 ). For two stars (HD24712 and HD 134214) the phase lag is 
consistent with the value of 0.25 at a 2-3<r level. 

How does one interpret a systematic quarter-period delay be- 
tween the variation of RV and line width? We recall that, in the 
absence of significant non-linearities and non-adiabatic effects in 
stellar pulsations, the velocity is shifted by 0.25 of the period with 
respect to radius changes. The negative extremum of (V) corre- 
sponds to the middle of the expansion part of the radius phase 
curve. At the maximum radius, RV returns to zero. Then the pos- 
itive part of the RV curve describes contraction phase, and ra- 
dius reaches minimum when (V) returns to zero after positive 
extremum. Adding to this picture the (V 2 ) variability discussed 
above, we find that line width changes approximately in antiphase 
(i.e., with a phase shift of 0.5) with respect to stellar radius varia- 
tion. In other words, the maximum contraction (minimum radius) 
phase point corresponds to maximum line width, whereas at the 
phase of maximum expansion (maximum radius) lines become nar- 
rower. This fundamental relation between pulsational modulation 
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Figure 9. Same as Fig. [8] for HD 122970, HD 128898, HD 137949 and 
HD 201601. 
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Figure 10. Illustration of the phase relation between pulsational modula- 
tion of stellar radius, radial velocity and line width discovered in our study. 
Symbols mark the phase point of maximum compression (filled circle) and 
maximum expansion (open circle). 
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Figure 11. Line width variation for the Nd III 6550 A line in the spectrum 
of 7 Equ. The average profile (dotted curve) is compared with the time- 
resolved spectra corresponding to minimum (solid curve) and maximum 
(dashed curve) of the second moment. 

of line width on the one hand and RV and radius changes on the 
other hand is sketched in Fig. 1101 

The tendency of REE lines to become broader at the maximum 
contraction phase is confirmed by investigation of the time-resolved 
spectra of 7 Equ. As was already mentioned above, our 2001 CFHT 
observations caught this star in a remarkably high-amplitude pulsa- 
tion state. The very high S/N and high resolution of these data re- 
veal pulsational profile variation synchronized with radius changes. 
Fig. II llcompares the average profile of the Nd III 6550 A line with 
the time-resolved spectra corresponding to nearly zero RV, but dif- 
ferent extrema of the {V 2 } phase curve. Time-resolved spectra for 
the minimum and maximum (V 2 ) exhibit different width, yet the 
equivalent width of the line does not change appreciably. This ev- 
idence confirms that our measurements of the second moment of 
REE in roAp stars diagnose pulsational perturbation of line width, 
not some other unrecognized line profile or line strength fluctua- 
tions. 



4.4 Short-term variation of moment amplitudes 

iKochukhov & Rvabchikova (2001b) discovered that the amplitude 
of the pulsational RV variability of the Nd III 6145 A line in a Cir 
(HD 128898) changes significantly within a time period of about 
one hour. This modulation of RV oscillations occurs on a time scale 
much shorter than the 4.47 1 d rotation period of a Cir and cannot 
be ascribed to beating between known close frequencies because 
photometric studies dKurtz et alJI 19941) show that this roAp star is 
essentially a monoperiodic pulsator. In a r ecent investigation of the 
pulsational RV variability of Prill lines. iKurtz et al. con- 
firmed the presence of short-term changes of RV amplitude in a Cir 
and demonstrated the existence of similar behaviour in a few other 
roAp stars. Among various possible explanations for the new spec- 
troscopic fr equencies, observe d in Prill but absent in broad-band 
photometry, 1 Kurtz et alj d2006h favoured the idea that at the forma- 
tion heights of Pr III lines the pulsational RV amplitude experiences 
stochastic variation due to growth and decay of the principal pulsa- 
tion mode. 

We have performed a new analysis of the short-term changes 
of pulsational variability using time series of both the first and the 
second moments measured in the UVES time-resolved spectra of 
a Cir. We explored the behaviour of several REE ions (Prill, Nd II, 



Ndll 5319 A 




o 2n 40 60 ao 100 

Time (mm) 
Tblll 5505 A 




-1.5 Li , , , i , , , i , , , i , , , i , , , i , , , i_: 

20 40 60 BO I00 120 

Tims (min) 



Figure 12. Modulation of RV amplitude for the Nd II 5319 A (upper panel) 
and Tb III 5505 A (lower panel) lines in HD 128898. Symbols show velocity 
measurements, solid line corresponds to the best-fitting superposition of the 
variability with Pi = 6.80 min and P2 = 7.34 min. The dashed curve 
illustrates temporal variation of the principal frequency amplitude. 



Table 6. Amplitude modulation for the oscillations in the first and second 
moments of REE lines in a Cir. The columns give line identification and the 
amplitude ratio V2/V1, which characterizes the strength of the secondary 
periodicity (P2 = 7.34 min) with respect to the amplitude of the principal 
frequency (Pi = 6.80 min). 



Ion A V2/V1 

(A) (V) (V 2 ) 



Pr III 


5284 


0.192 ±0.018 


0.331 ±0.026 


Ndm 


5286 


0.154 ±0.025 


0.392 ±0.062 


Ndm 


5294 


0.255 ±0.031 


0.360 ± 0.046 


Prill 


5300 


0.165 ±0.022 


0.369 ±0.029 


Ndll 


5319 


0.073 ± 0.022 


0.224 ± 0.053 


Tblll 


5505 


0.220 ±0.018 


0.371 ±0.023 


Ndm 


5802 


0.141 ±0.024 


0.423 ± 0.059 


Tblll 


5847 


0.188 ±0.038 


0.334 ± 0.050 


Ndm 


5851 


0.198 ±0.024 


0.370 ±0.054 



Nd III, Tb III) and considered the results in the light of the newly 
discovered depth dependence of the pulsational LPV pattern and 
line width oscillations. Analysis of the amplitude variability for 
the REE line profile moments was based on the least-squares fit- 
ting introduced in Sect. 14.31 except that here we adjust amplitude 
and phase for the principal frequen cy (Pi = 6.80 mi n) and a sec- 
ondary periodicity P2 = 7.34 min dKurtz et alj |2006). responsible 
for the observed short-term modulation. The strength of the ampli- 
tude modulation is then characterized by the respective amplitude 
ratio V2/V1, estimated separately for the first and second moments. 

Numerical results of the least-squares analysis of the moment 
time series is summarized in Table [6] for nine REE lines. Clearly, 
we find different amounts of RV amplitude modulation for differ- 
ent REE ions. The secondary frequency is at the detection threshold 
for the Nd II 5319 A line. The V2/V1 amplitude ratio increases by a 
factor of 2-3 towards the higher atmospheric layers where doubly 
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Table 3. Results of the least-squares fitting of the first and second moment variation of the REE lines in HD9289, HD 12932, HD 19918, HD 24712 and 
HD 101065. The columns give line identification, amplitude and phase (Vu, ipu) for the (V) variation at the main frequency and amplitude of the (V) 
variability (V12) at the first harmonic. The same quantities are listed for the second moment (V 2 ). Amplitudes are measured in kms -1 for (V) and in 
km 2 s -2 for (V 2 ), phases are given in units of the pulsation period. The HD number of the star is followed by the phase difference (fi2l — 'fin, estimated 
using the lines of Pr III and Tb III. 



Ion A (V) (V 2 ) 

(A) Vndcms- 1 ) ip n Viatkms- 1 ) V 2 i (km 2 s~ 2 ) ^21 V22 (km 2 s" 2 ) 
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Table 4. Results of the least-squares fitting of the first and second moment variation of the REE lines in HD 122970, HD 128898, HD 134214 and HD 137949. 
See caption of Table|5]for the explanation of this table's content. 



Ion A (V) (V 2 ) 

(A) Vnikms- 1 ) tp n V12 (kmr 1 ) V 2 i (km 2 s~ 2 ) <p 2 l V22 (km 2 s~ 2 ) 



















HD 122970, 


(ip 2 l - ipll) = 


0.22 ± 0.04 
























Prill 


5284 





.614 


± 


0. 


016 


0. 


.891 ± 0.004 


0.032 ± 0.016 


0.930 ±0. 


.068 





.926 


± 


0. 


012 


0. 


465 


± 


0. 


.068 


Ndlll 


5286 





.535 


± 


0. 


.027 


0. 


.844 ± 0.008 


0.021 ±0.027 


0.309 ± 


.120 





.729 


± 


0. 


063 


0. 


.227 


± 


0. 


.121 


Ndm 


5294 





.434 


± 


0. 


011 





.853 ± 0.004 


0.025 ±0.011 


0.270 ±0. 


.040 





.875 


± 


0. 


024 


0. 


118 


± 





.040 


Prill 


5300 





.586 


± 


0. 


.014 


0. 


.886 ± 0.004 


0.020 ±0.014 


0.283 ± 0. 


.058 





.086 


± 


0. 


033 


0. 


245 


± 


0. 


.058 


Ndll 


5311 





.490 


± 


0. 


023 


0. 


.842 ± 0.007 


0.028 ± 0.023 


0.103 ±0. 


.101 


0. 


.828 


± 


0. 


154 


0. 


154 


± 


0. 


.101 


Ndll 


5319 





.430 


± 


0. 


012 





.843 ± 0.005 


0.006 ± 0.012 


0.074 ±0. 


.041 





.206 


± 


0. 


089 


0. 


163 


± 





.041 


Tbm 


5505 


1 


.107 


± 


0. 


026 


0. 


.953 ± 0.004 


0.136 ±0.026 


1.504 ±0. 


.146 





.213 


± 


0. 


.015 


0. 


838 


± 





.145 


Ndm 


5802 





.511 


± 


0. 


020 





.846 ± 0.006 


0.019 ±0.020 


0.126 ±0. 


.085 





.953 


± 


0. 


.108 


0. 


143 


± 





.085 


Tbm 


5847 


1. 


.095 


± 


0. 


039 


0. 


.930 ± 0.006 


0.113 ±0.038 


1.383 ±0. 


.210 


0. 


.113 


± 


0. 


024 


1. 


.185 


± 


0. 


.211 


Ndm 


5851 





.515 


± 


0. 


019 


0. 


.853 ± 0.006 


0.030 ±0.019 


0.553 ±0. 


.096 


0. 


.948 


± 


0. 


028 


0. 


422 


± 


0. 


.096 


Tbm 


6323 


1 


.197 


± 


0. 


053 





.946 ± 0.007 
HD 128898, 


0.169 ±0.052 

(ip 2 i - (fin) = 


1.239 ±0. 
0.28 ± 0.04 


.320 





.176 


± 


0. 


.041 


1. 


096 


± 





.315 


Prill 


5284 





.509 


± 


0. 


011 


0. 


.670 ± 0.003 


0.017 ±0.011 


2.715 ±0. 


.084 





.887 


± 


0. 


005 


0. 


046 


± 


0. 


.084 


Ndm 


5286 





.708 


± 


0. 


.018 


0. 


.359 ± 0.004 


0.075 ±0.018 


2.248 ± 0. 


.132 





.687 


± 


0. 


009 


0. 


346 


± 





.132 


Ndm 


5294 





.287 


± 


0. 


009 





.604 ± 0.005 


0.025 ±0.009 


1.740 ±0. 


.080 





.911 


± 


0. 


.007 


0. 


002 


± 





.080 


Prill 


5300 





.513 


± 


0. 


012 


0. 


.660 ± 0.004 


0.020 ±0.012 


2.916 ±0. 


.100 





.937 


± 


0. 


005 


0. 


133 


± 


0. 


.099 


Ndll 


5319 





.815 


± 


0. 


.018 


0. 


.333 ± 0.003 


0.067 ±0.018 


2.553 ± 0. 


.132 


0. 


.778 


± 


0. 


008 


0. 


289 


± 


0. 


.131 


Tbm 


5505 





.777 


± 


0. 


017 


0. 


.840 ± 0.003 


0.059 ±0.017 


5.088 ± 0. 


.153 





.137 


± 


0. 


005 


0. 


014 


± 


0. 


.154 


Ndm 


5802 





.512 


± 


0. 


013 


0. 


.432 ± 0.004 


0.036 ± 0.013 


1.930 ± 0. 


.110 





.798 


± 


0. 


009 


0. 


236 


± 





.110 


Tbm 


5847 





.854 


± 


0. 


033 





.834 ± 0.006 


0.076 ± 0.033 


5.376 ± 0. 


.257 





.123 


± 


0. 


008 


0. 


295 


± 


0. 


.257 


Ndm 


5851 





.563 


± 


0. 


014 


0. 


.441 ± 0.004 


0.046 ± 0.014 


2.088 ± 0. 


.105 





.769 


± 


0. 


008 


0. 


431 


± 


0. 


.105 


Tbm 


6323 





.746 


± 


0. 


.040 


0. 


.802 ± 0.008 
HD 134214, 


0.055 ±0.040 

(ip 2 l - <pil) = 


3.781 ±0. 
0.20 ±0.03 


.271 





.108 


± 


0. 


011 


0. 


612 


± 


0. 


.270 


Prill 


5284 





.387 


± 


0. 


016 


0. 


.226 ± 0.007 


0.010 ±0.016 


2.241 ±0. 


.102 





.376 


± 


0. 


007 


0. 


063 


± 


0. 


.103 


Ndm 


5286 





.553 


± 


0. 


.021 


0. 


.921 ± 0.006 


0.014 ±0.021 


1.233 ± 0. 


.127 


0. 


.099 


± 


0. 


016 


0. 


100 


± 


0. 


.128 


Ndm 


5294 





.207 


± 


0. 


.010 


0. 


.095 ± 0.007 


0.010 ±0.010 


0.728 ± 0. 


.052 


0. 


.374 


± 


0. 


011 


0. 


052 


± 


0. 


.052 


Prm 


5300 





.362 


± 


0. 


015 





.205 ± 0.006 


0.013 ±0.015 


2.266 ± 0. 


.104 





.425 


± 


0. 


007 


0. 


122 


± 





.104 


Ndn 


5311 





.508 


± 


0. 


019 


0. 


.894 ± 0.006 


0.001 ±0.020 


0.482 ± 0. 


.073 





.120 


± 


0. 


024 


0. 


.242 


± 


0. 


.072 


Ndll 


5319 





.522 


± 


0. 


.018 


0. 


.928 ± 0.006 


0.010 ±0.018 


0.843 ± 0. 


.053 





.214 


± 


0. 


010 


0. 


155 


± 


0. 


.053 


Tb III 


5505 





.508 


± 


0. 


027 


0. 


.410 ± 0.009 


0.026 ± 0.027 


2.713 ± 0. 


.159 





.614 


± 


0. 


009 


0. 


143 


± 


0. 


.158 


Ndm 


5802 





.444 


± 


0. 


.018 


0. 


.963 ± 0.006 


0.004 ±0.018 


1.137 ±0. 


.085 


0. 


.155 


± 


0. 


012 


0. 


137 


± 


0. 


.085 


Tbm 


5847 





.549 


± 


0. 


043 





.395 ± 0.013 


0.017 ±0.043 


4.243 ± 0. 


.375 





.613 


± 


0. 


.014 


0. 


513 


± 





.377 


Ndm 


5851 





.490 


± 


0. 


019 





.964 ± 0.006 


0.019 ±0.019 


1.308 ± 0. 


.098 





.179 


± 


0. 


012 


0. 


116 


± 





.099 


Tbm 


6323 





.673 


± 


0. 


060 


0. 


.413 ±0.014 
HD 137949, 


0.045 ± 0.059 

{ip 2 l - ifn) = 


3.683 ±0. 
0.27 ±0.07 


.546 


0. 


.606 


± 


0. 


023 


0. 


527 


± 


0. 


.549 


Prill 


5284 





.056 


± 


0. 


.004 


0. 


.307 ±0.010 


0.028 ± 0.004 


0.535 ±0. 


.027 





.500 


± 


0. 


008 


0. 


.242 


± 





.027 


Ndm 


5286 





.023 


± 


0. 


005 





.723 ± 0.036 


0.059 ± 0.005 


0.383 ± 


.060 





.246 


± 


0. 


025 


0. 


588 


± 





.059 


Ndm 


5294 


0. 


.128 


± 


0. 


006 


0. 


.336 ± 0.007 


0.011 ±0.006 


0.779 ± 0. 


.064 


0. 


.582 


± 


0. 


013 


0. 


246 


± 


0. 


.064 


Prm 


5300 





.069 


± 


0. 


005 





.253 ± 0.011 


0.035 ±0.005 


0.576 ± 0. 


.044 





.492 


± 


0. 


.012 


0. 


293 


± 


0. 


.044 


Ndn 


5311 





.096 


± 


0. 


005 





.715 ±0.007 


0.047 ± 0.005 


0.135 ±0. 


.020 





.966 


± 


0. 


024 


0. 


499 


± 





.020 


Ndn 


5319 





.114 


± 


0. 


005 


0. 


.795 ± 0.008 


0.066 ± 0.005 


0.157 ±0. 


.034 


0. 


.184 


± 


0. 


.035 


0. 


642 


± 


0. 


.034 


Ndm 


5802 





.064 


± 


0. 


006 


0. 


.298 ± 0.016 


0.042 ± 0.006 


0.585 ± 0. 


.065 





.330 


± 


0. 


.018 


0. 


.228 


± 


0. 


.065 


Tbm 


5847 


0. 


.096 


± 


0. 


009 


0. 


.003 ± 0.015 


0.062 ± 0.009 


0.268 ± 0. 


.072 





.355 


± 


0. 


.043 


0. 


587 


± 


0. 


.072 


Ndm 


5851 





.064 


± 


0. 


008 





.244 ± 0.020 


0.070 ± 0.008 


0.604 ± 0. 


.103 





.443 


± 


0. 


.027 


0. 


.532 


± 





.103 


Tbm 


6323 





.175 


± 


0. 


017 





.051 ± 0.016 


0.099 ±0.017 


0.971 ±0. 


.169 





.357 


± 


0. 


028 


0. 


929 


± 





.169 



ionized REE lines are formed. This is further illustrated in Fie. 1121 
where we show the observed and best-fitting RV curve for the Nd II 
5319 A and Tbm 5505 A lines. One can clearly see that pulsa- 
tions in the latter line change appreciably with time, whereas much 
weaker modulation is observed for the former line. This figure and 
Table[6]provides the first definite demonstration of the height vari- 
ation of the short-term amplitude modulation effect in a Cir. More- 
over, time-series analysis of the second moment reveals a much 
stronger (by a factor 1.5-3) amplitude modulation compared to the 
RV results. This is unexpected because the relative amplitudes of 
the two frequencies should be the same for both (V) and (V 2 ) 
( Math ias et all 19941) . 



Thus, it is tempting to speculate that, whatever effect gives 
rise to the observed RV amplitude variability in the upper atmo- 
sphere of a Cir, it acts more strongly on the line width variation 
than on the first moment. Considering these results in the light of 
our discovery of the steep growth of the line width oscillations with 
height, we suggest that the observed short-term modulation is actu- 
ally a property of the physical mechanism responsible for the line 
width and asymmetric LPV behaviour, while the short-term am- 
plitude changes visible in RV represent a secondary consequence 
of the action of this hypothetical mechanism. Further investigation 
of several roAp stars is required to verify this hypothesis. In par- 
ticular, a full night of uninterrupted time-resolved monitoring of 
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Table 5. Results of the least-squares fitting of the first and second moment variation of the REE lines in HD 201601 . See caption of Table[5]for the explanation 
of this table's content. The phase difference <£>2i — ¥>n is obtained using the lines of Pr III, Nd III and Tb III. 



Ion A (V) (V 2 ) 

(A) Vndcms" 1 ) ip u V^fkms" 1 ) V 2 i (km 2 s" 2 ) V2 i V22 (km 2 s~ 2 ) 





















HD201601 1 , (ip2i - ¥>il 


} = 0.21 






















Ndlll 


6145 





.458 


± 


0. 


016 





.178 ± 


0.005 


0.056 ±0.016 


2.989 ±0.186 


0. 


.390 


± 


0. 


010 


0. 


574 


± 





.194 


Prill 


6160 





.642 


± 


0. 


.018 





.179 ± 


0.004 


0.067 ±0.019 


3.192 ±0.183 





.385 


± 


0. 


009 


0. 


302 


± 





.188 


















HD 


201601 2 . 


, {ip2l - ¥>ll) = 


0.26 ±0.05 






















Ndlll 


6550 


0. 


.737 


± 


0. 


011 


0. 


.605 ± 


0.002 


0.060 ±0.011 


5.038 ±0.248 





.836 


± 


0. 


008 


0. 


066 


± 


0. 


.246 


Ndll 


6650 


1 


.131 


± 


0. 


.021 


0. 


.399 ± 


0.003 


0.130 ±0.020 


2.063 ±0.178 


0. 


.607 


± 


0. 


.014 


0. 


969 


± 


0. 


.177 


Ndlll 


6691 


1 


.078 


± 


0. 


020 





.509 ± 


0.003 


0.097 ±0.020 


2.792 ±0.280 


0. 


.791 


± 


0. 


016 


1. 


193 


± 





.276 


Prill 


6706 





.841 


± 


0. 


.027 





.599 ± 


0.005 


0.070 ± 0.027 


4.495 ±0.392 


0. 


.901 


± 


0. 


.014 


1. 


.081 


± 





.393 


















HD 


201601 3 . 


, (tp2l — t<9ll) = 


0.26 ±0.06 






















Prill 


5284 


0. 


.158 


± 


0. 


006 


0. 


.345 ± 


0.006 


0.007 ± 0.006 


0.405 ± 0.064 





.581 


± 


0. 


.025 


0. 


075 


± 


0. 


.064 


Ndlll 


5286 





.319 


± 


0. 


013 





.185 ± 


0.007 


0.011 ±0.013 


0.714 ± 0.144 





.431 


± 


0. 


032 


0. 


079 


± 





.143 


Ndlll 


5294 





.104 


± 


0. 


007 





.372 ± 


0.011 


0.002 ± 0.007 


0.791 ±0.067 


0. 


.564 


± 


0. 


.014 


0. 


101 


± 





.068 


Prill 


5300 


0. 


.234 


± 


0. 


.010 


0. 


.346 ± 


0.006 


0.025 ±0.010 


1.094 ±0.114 





.577 


± 


0. 


017 


0. 


152 


± 


0. 


.116 


Ndll 


5311 


0. 


.277 


± 


0. 


013 


0. 


.139 ± 


0.008 


0.032 ±0.013 


0.253 ±0.056 





.473 


± 


0. 


036 


0. 


085 


± 


0. 


.056 


Ndll 


5319 





.310 


± 


0. 


.010 





.178 ± 


0.005 


0.018 ±0.010 


0.303 ± 0.042 


0. 


.465 


± 


0. 


022 


0. 


095 


± 





.042 


TbiII 


5847 





.222 


± 


0. 


038 





.499 ± 


0.028 


0.071 ±0.038 


1.193 ±0.349 


0. 


.851 


± 


0. 


047 


0. 


268 


± 





.353 


Ndlll 


5851 


0. 


.218 


± 


0. 


.016 


0. 


.375 ± 


0.012 


0.028 ± 0.016 


0.667 ±0.130 





.661 


± 


0. 


.031 


0. 


290 


± 


0. 


.129 




















HD201601 4 , <v3 2 i - V11 


> = 0.15 






















TbiII 


5847 


0. 


.178 


± 


0. 


.046 


0. 


.823 ± 


0.041 


0.034 ± 0.046 


1.318 ±0.489 





.997 


± 


0. 


059 


0. 


.440 


± 


0. 


.481 


Ndlll 


5851 





.199 


± 


0. 


016 





.705 ± 


0.013 


0.012 ±0.016 


0.515 ±0.144 


0. 


.838 


± 


0. 


044 


0. 


.250 


± 





.144 




















HD201601 5 , (<p 2 i - Vii 


> = 0.25 






















Ndlll 


6550 


0. 


.118 


± 


0. 


013 


0. 


.906 ± 


0.017 


0.002 ± 0.013 


0.579 ±0.167 





.182 


± 


0. 


046 


0. 


231 


± 


0. 


.165 


Ndll 


6650 





.149 


± 


0. 


028 





.715 ± 


0.031 


0.015 ± 0.029 


0.214 ±0.166 


0. 


.933 


± 


0. 


123 


0. 


173 


± 





.166 




















HD201601 6 , (<p 2 i - Vii 


> = 0.23 






















Ndlll 


6145 


0. 


.155 


± 


0. 


008 


0. 


.521 ± 


0.008 


0.013 ± 0.008 


0.843 ±0.148 





.751 


± 


0. 


028 


0. 


102 


± 


0. 


.147 


Prill 


6160 


0. 


.175 


± 


0. 


011 


0. 


.546 ± 


0.010 


0.025 ±0.011 


0.910 ±0.082 





.783 


± 


0. 


014 


0. 


.188 


± 


0. 


.081 



a Cir would be extremely useful to discriminate between a truly 
stochastic and a quasiperiodic nature of the variation of oscillation 
amplitudes in the upper atmospheric layers. 



5 INTERPRETATION 

Our survey of the pulsational line profile variability in sharp-lined 
roAp stars has revealed a new type of LPV behaviour and has 
demonstrated its significant depth dependence. In summary, our ob- 
servational findings are: 

(i) The asymmetric, blue-to-red wave pattern is common in 
time-resolved residual spectra of roAp stars. This behaviour is more 
pronounced for REE lines formed in the higher atmospheric layers, 
whereas lines originating at the bottom of the REE-enriched cloud 
often display S-shaped, symmetric, blue-to-red-to-blue LPV pat- 
tern. 

(ii) The doubly ionized REE lines from the uppermost part of 
the stellar atmosphere show strong single-wave pulsational vari- 
ability of line width. These oscillations typically occur in antiphase 
with the pulsational changes of stellar radius. They also correlate 
with the presence of asymmetric LPV and exhibit the same height 
dependence. 

(iii) The efficiency of the mechanism responsible for the asym- 
metric LPV pattern and single-wave line width variation in 
HD 128898 changes on relatively short time scales. 

These results, in particular the new type of pulsational LPV 
described by (i) and (ii), challenge our understanding of pulsations 
in roAp stars and call for theoretical explanation. In this section 



we present a series of line profile calculations for oblique non- 
radial pulsations and compare observations with theoretical expec- 
tations. We show that the standard models including only spherical 
harmonic perturbations of velocity fail to reproduce the observed 
line profile and moment variability. This led us to introduce a new 
oblique pulsator model, which accounts for additional line width 
oscillations. Using this new modelling framework, we are able to 
explain all the main characteristics of the pulsational LPV in sharp- 
lined roAp stars. 



5.1 Spectrum synthesis calculations for oblique pulsators 

Periodic variation of velocity in non-radial pulsators is described 
with the spherical harmonics Y™ (8, (f>) and their derivatives: 



(3) 



KV V 



-YP 



where u is the angular pulsation frequency, V p is the pulsation am- 
plitude and K is the ratio of the vertical to horizontal amplitudes. 
Non-radial oscillations in roAp stars are oblique. Therefore, spher- 
ical angular coordinates 8, (f> refer to a reference system whose ori- 
entation is defined by the angles /3 and \ with respect to the stellar 
axis of rotation. In the classical oblique pulsator model, the pulsa- 
tion axis is assumed to be aligned with the axis of a quasi-dipolar 
magnetic field. In this case, 8 and (f> are spherical angular coordi- 
nates in the magnetic reference frame. 

Eq. l[3} can easily be generalized (see lKochukhovll2004al) for 
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the situation when pulsation modes are distorted by stellar rotation 
and/or a strong magnetic field and are described by superposition 
of several spherical harmonic functions. 

The horizontal to vertical amplitude ratio is usually approxi- 
mated as 



K = 



GM 
17W 



(4) 



(e.g.. ISmevers & Tassoull ["l987t) . This relation predicts that ver- 
tical oscillations should be dominant (i.e., K <C 1) for high- 
overtone p-modes in roAp stars. Consequently, horizontal displace- 
ment was often neglected in previous studies of roAp pulsations 
(Bigot & Dziembowski 2002). However, the distorted magnetoa- 
coustic pulsations modes in the theory of ISaio & Gautschvl (2004) 
show significant horizontal amplitudes (K ~ 1) in the lower part 
of stellar atmosphere. 

For the purpose of computing theoretical LPV, the pulsation 
velocity field described by Eq. ^ is transformed to the Descartes 
coordinate system of the observer. The line of sight velocity com- 
ponent, V°, which enters spectrum synthesis, combines the contri- 
bution of pulsations and the Doppler shift due to solid body stellar 
rotation with the projected rotational velocity v c sin i. For the sake 
of brevity here we omit details of coordinate transformations, def- 
inition of spherical harmonics and calculation of their der i vatives . 
Interested readers can find this information in iKochukhovl d2004j. 
l2005t) . 

Even neglecting possible rotational and magnetic distortion of 
the non-radial pulsations, the model formulated above has a large 
number of free parameters. One has to specify the amplitude and 
structure of the pulsation mode, V p ,K,i, m, as well as the orienta- 
tion of the pulsation axis f3, \, and the stellar v c sin i and inclination 
angle i. In the limit of negligible v a sini applicable to many roAp 
stars, the angles /3, x and * can be replaced by a single parameter a 
- the angle between the line of sight and pulsation axis. 

In order to be able to explore properties of LPV for many dif- 
ferent pulsation geometries in a reasonable amount of time, we ap- 
proximate the local line profile with a Gaussian function, charac- 
terized by a constant central depth D and full width at half maxi- 
mum W. Then the disk-integrated line profile S is obtained by the 
weighted summation over all visible surface zones 



D 



* -41n2 

x r, e w 



where the weight function 

3(1 -u + Ufx^HjAj 



Ti 



7r(3 — u) 



(5) 



(6) 



accounts for unequal projected areas of surface elements and in- 
cludes a linear limb-darkening law with coefficient u, and /x de- 
notes the cosine of the angle between the line of sight and the sur- 
face normal. In all simulations described below we use D = 0.5, 
W = 10 kms" 1 , u — 0.5. V p and K are chosen in such a way 
that the resultant RV amplitude is compatible with the one typi- 
cally found in roAp stars (~ 1 kms -1 ). The stellar surface is sam- 
pled on a 3909-element g rid with unequal number of su rface zones 
at different latitudes (see lPiskunov & Kochukhovll2002h . Calcula- 
tions are performed for a set of 20 equally-spaced pulsation phases. 
Eq. |QJ is used to infer variation of the first and second moments 
from theoretical line profile time series. 

Our LPV synthesis model is very similar to the one usually 
used in a broader context of spectral variability studies of non- 



radial pulsations in sl owly rotating stars dAerts & Waelkensll 19931 : 
ISchriivers etafll 19971) . Using this standard description of oscilla- 
tions and accounting for the additional complication caused by 
oblique mode geometry, we carried out an extensive set of com- 
putations aiming to find a combination of input parameters that 
produces LPV compatible with observations. Pulsational LPV and 
moment variability was investigated for 

(1) all possible modes with I ^ 3; 

(2) projected rotational velocities in the range 0-10 km s _1 ; 

(3) modes with significant (K ~ 1) and dominant (K ^> 1) 
horizontal fluctuations; 

(4) distorted dipolar modes give n by superposition of £ = 1 , 
m = —1,0, 1 spherical harmonics teigot & Dziembowskj||20o"2h 
or by superposition of £ = 1 and £ — 3 axisymmetric harmonics 
dSaio & Gautschvll2004l) : 

(5) modes containing two pulsation components with different 
K and a phase shift relative to each other - this model is meant 
to simulate superposition of magnetic and acoustic running wave 
comp onents possi bly existing in the upper atmospheres of roAp 
stars JCunhiJl2006l) . 

We found that none of the pulsation structures studied leads to 
a qualitative agreement with observations. A typical failure of the 
theoretical model is illustrated in Fig. 1 13b for the I = 1, m = 
mode viewed from the pulsation pole. The residual profile varia- 
tion has a symmetric S-shape and the line width changes have low- 
amplitude, double-wave character. In the absence of stellar rota- 
tion, LPV for all pulsation modes have these basic properties - the 
same conclusion follo ws from the theoretical analysis of the second 
moment behaviour ( Aerts et al. 1992; Kochukhov 2005). This vari- 
ability is not incompatible with observations of weak Nd II lines in 
a few roAp stars, but clearly contradicts pulsation signatures found 
for the majority of doubly ionized REE lines. 

When significant (v e sini ^> V p ) rotational velocity is intro- 
duced in the model and the pulsation pole of the dipolar mode is 
offset from the line of sight, one may encounter LPV character- 
ized by significant single-wave line width changes. However, the 
phase shift of RV and line width oscillations is inconsistent with 
observations. Moreover, for this pulsational geometry a small RV 
amplitude corresponds to a quite strong LPV, which has a fully 
symmetric (both in time and relative to the line centre) H-shape 
and thus appears to be very different from the observed blue-to-red 
residual profile shift. 

Our calculations show that all types of complex modes, ob- 
tained by combining different spherical harmonics or by combining 
components with different phase shifts, also fail to achieve a satis- 
factory agreement with observations. The key features of the roAp 
line profile behaviour - single-wave line width oscillations with a 
characteristic phase shift and asymmetric LPV pattern - are never 
found for the theoretical models explored. Thus, we conclude that 
the standard picture of the velocity field in non-radial pulsators and 
trivial extensions of this picture are inadequate for the REE lines 
which form in the upper atmospheres and show the most promi- 
nent pulsation signatures in roAp stars. 

Admittedly, the simplified line formation model adopted in 
our calculations does not allow one to properly explore effects of 
the depth dependent pulsation amplitude and phase. However, we 
expect that any superposition of several harmonic oscillations will 
give another harmonic oscillation. Consequently, variation of the 
pulsation wave properties with height will not introduce any qual- 
itatively new behaviour and will not help to explain the puzzling 
spectroscopic observations of roAp stars. 
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Figure 13. Line profile variation of a non-rotating non-radial pulsator. a) Spectrum variability for the i = 1, m = 0, V p = 2.0 km s~ x , K = mode viewed 
from the pulsation pole. The local profile is represented by a constant Gaussian with FWHM W = 10 km s . b) Effect of adding harmonic variability of the 
line width with an amplitude 8W = 1.5 kms -1 and a phase shift of 0.25 with respect to the pulsational velocity variation. In both cases the star is assumed 
to have v e sin i = km s ~~ 1 . In each panel the left plot shows the average line profile on top and time series of the difference spectra (covering two pulsation 
cycles) below. The right panels illustrates variation of the first (upper plot) and second (lower plot) line profile moments. 
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Figure 14. The same as in Fig.[6p for the oblique non-radial pulsator with a) v e sini = 10 kms 1 and b) it sini = 20 kms 1 . 



It appears that the key problem lies in the formulation of the 
standard pulsation framework: the time dependent factor in Eq. l[3) 
postulates symmetry of the compression and expansion parts of 
pulsation cycle. In reality the opposite, asymmetric behaviour is 
observed in roAp stars (see Sect.|4]and Fig. Hot. 



5.2 A new oblique pulsator model 

We have discovered that the enigma of LPV in roAp stars can be 
resolved with a simple modification of the classical oblique pul- 
sator framework. Instead of looking for a combination of pulsation 
parameters which gives the variation of the disk-integrated quan- 
tities illustrated in Fig. [10] we suggest that this empirical relation 
between radius, RV and second moment changes actually reflects 
a local periodic variation of velocity and line width. Thus, the key 
idea of our new oblique pulsator model is that the pulsational mech- 
anism responsible for the RV variability of REE lines gives rise to 
additional, sinusoidal changes of line width. This variability occurs 
approximately in antiphase with the pulsational radius changes and 
hence shows a quarter of period phase shift with respect to RV os- 



cillations. Then the LPV expected from this model is a combination 
of two effects: 

(i) the usual pulsational velocity changes described by spherical 
harmonic perturbations and 

(ii) sinusoidal variability of line profile width, synchronized 
with changes of stellar radius. 

It is reasonable to assume that the surface distribution of the 
line width amplitude is similar to the vertical pulsational displace- 
ment and is described by 

w(t,e,(f>) ^Wo + swYne,^' 1 ^-^. a) 

The theoretical LPV can be calculated by incorporating this time 
and position-dependent line width in the disk-integration formula 
(|5). A striking effect of introducing variable line width in the line 
profile synthesis is demonstrated by Fig.|13b. The pulsation model 
used for this figure is identical to the model of Fig. 113b . except 
that now line width changes with an amplitude 8W = 1.5 kms - 
around Wo = 10 kms -1 and a phase shift ipw = 0.25. This 
leads to a remarkable agreement between theoretical residual pro- 
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Figure 15. Dependence of the line profile variation pattern on a) amplitude 8W of the line width changes, b) phase shift Aip between the line width and 
pulsation velocity variability, and c) the angle a between the pulsation axis and the line of sight. The standard set of the oblique pulsator parameters is i = 1, 
m = 0, V p = 2.0 kms" 1 , K = 0, SW = 1.5 kins -1 , a = 0°. In each subpanel the average line profile and the standard deviation (arbitrary scale) are 
presented on top. The time-resolved difference spectra, covering two pulsation cycles, are plotted below. 



file variation and observations of roAp stars. Model spectra show 
smooth, blue-to-red moving pattern, which is very similar to the 
LPV discussed in Sect. 14.21 The additional variation of line width 
does not distort the sinusoidal behaviour of the RV. Considering the 
second moment oscillations, we find that a large-amplitude single- 
wave variability superimposed on top of a small-amplitude double- 
wave curve results in nearly sinusoidal variation of ( V 2 } . The phase 
shift assumed for the time dependence of the local velocity and line 
width propagates to the disk-integrated observables, allowing us to 
reproduce the observed phase lag of (V 2 ) with respect to (V). 



An interesting modification of the LPV pattern calculated with 
our model is associated with increasing v e sin i. The LPV corre- 
sponding to u e sini = kms -1 is presented in Fig. 113b. Fig. 1 14b 
and b show LPV for the same oblique pulsator geometry and ampli- 
tude, but with v e sin i increased to 10 and 20 kms -1 , respectively. 
The blue-to-red shift is still discernible for the v e sin i = 10 km s _1 
model, although it is not as smooth as in nonrotating star. In the 
case of rapid rotation (Fig. 114b). the signature of line width vari- 
ability disappears and the LPV shows an almost symmetric pat- 
tern, as one would expect for the standard oblique pulsator. Thus, 
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if the SW amplitude adopted for these simulations is realistic, we 
predict that a signature of the pulsational variation of line width 
will be visible only in slowly rotating stars (v e sin i < 10 km s -1 ), 
but not in relatively rapid rotators (v e sini>20 kms -1 ). In fact, 
this prediction is verified by observations. We find asymmetric 
LPV in all slowly rotating roAp stars, including HD9289 and 
HD 128898 (u c sini~10 fans" 1 ), whereas the study of high- 
amplitude, rapidly rotating (u c sini = 33 kms -1 ) roAp pulsator 
HD 83368 reveals a symmetric LPV pattern and no evidence of 
blue- to-red shifting features in residual line profiles (Kochukhov 
l2006h . 

We have performed additional calculations with the aim of ex- 
amining the sensitivity of the theoretical LPV to the parameters 
of the new oblique pulsator model. Fig. Q3] illustrates the modifi- 
cations of profile variations due to changes of the amplitude and 
phase of the line width oscillations and due to changes of the an- 
gle a between the line of sight and pulsation axis. Gradual increase 
of the line width amplitude transforms symmetric, S-shaped LPV 
into asymmetric, blue-to-red pattern (see Fig.|15b). Thus, the diver- 
sity and height dependence of the LPV observed in roAp stars can 
be plausibly attributed to different V P -AW combinations and to a 
growth of the line width amplitude with height. Qualitative compar- 
ison of the simulated profile time series and observations suggests 
that for most stars SW lies in the range from to 1.5-2 kms -1 . 
Fig. 1 1 5 fa also demonstrates that the shape of the standard devia- 
tion profile changes from a double peak to single peak when SW is 
increased. This agrees with our observations of the correlation be- 
tween LPV pattern and the standard deviation shape. The fact that 
we never observe a dominant central peak in the standard deviation 
signatures also points to SW amplitudes below 1.5-2 kms -1 . 

Changing the ipw parameter from to ~ 0.15 also strongly 
affects the LPV character. Fig. |15b clearly shows that it is a non- 
zero phase shift between the velocity and line width oscillation that 
introduces an asymmetry in LPV. The blue-to-red shifting features 
are obtained for positive <pw exceeding ~ 0.1. On the other hand, 
negative ipw results in the opposite, red-to-blue, pattern, possibly 
detected in Nd II lines of HD 9289. The morphology of the LPV is 
relatively unaffected by variation of the phase shifts between 0.15 
and 0.35. Interestingly, a ipw different from and 0.25 produces 
an asymmetric standard deviation, which is sometimes shifted with 
respect to the line centre. This provides the first explanation of the 
origin of distorted and shifted standard deviation profiles frequently 
observed in slowly rotating roAp stars. 

Rotational modulation of the LPV pattern presented in 
Fig. |15t shows a small change of profile variation as the line of 
sight shifts from the pulsation pole to the pulsation equator of the 
oblique dipolar mode. A dramatic change of the LPV is associated 
with a narrow range of rotation phases around a — 90°, when 
the mode is viewed equator-on. For this pulsation geometry the RV 
variation is negligible and LPV is caused entirely by the line width 
oscillations. 

The line profile variation calculations presented in this sec- 
tion convincingly demonstrate the success of our new oblique pul- 
sator model. Simultaneous pulsational variation of velocity and line 
width naturally explains a wide range of the LPV behaviour ob- 
served in roAp stars. We showed that, by adjusting a few parame- 
ters (amplitude and phase of the line width variation) of our model, 
it becomes possible to reproduce the diverse variations found for 
different groups of REE lines. 
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Figure 16. Line profile variation expected for the oblique non-radial 
pulsator within the fra mework of the shock wave model proposed by 
IShibahashi et alj fe004l) . The velocity field geometry is that of the pole-on 
i = 1, m = 0, Vp = 1.6 kms -1 pulsation mode. The temporal varia- 
tion is given by the smoothed saw-tooth function n_1 ex P (' na; *)- 
Other parameters and the layout of the figure are the same as in Fig.EIJi. 

5.3 A shock wave signatures in the LPV of roAp stars? 

Recently, another interpretation of pulsational pr o file v ariations 
in roAp stars was proposed by IShibahashi etall d2004l) . These 
aut hors suggested that the smooth, blue-to-red shift observed 
by iKochukhov & Rvabchikovj d2001ah in the LPV of the Nd III 
6145 A line in 7 Equ represents the si gnature of shock w aves prop- 
agating in the high atmospheric layers. [shibahashi"et al. deduce the 
velocity of the pulsation wave, not from the amplitude of the RV 
oscillation (~ 500 m s _1 for this Nd III line), but from the full span 
of the residual profile variation (see Fig.QJ. This gives a maximum 
speed of ~ 18 km s , which is s ubstantially higher than the sound 
speed. Based on these arguments, Shibaha shi et alj claim that LPV 
is affected by shocks and the harmonic time dependence in Eq. l[3} 
should be replaced by a saw-tooth law typical for large-amplitude 
non-linear stellar oscillations. 

W e believe that the basic assumption of Shib ahashi et al.1 
d2004[) is incorrect because it is based on an erroneous interpre- 
tation of the pulsational profile distortions in terms of the physical 
velocity of the pulsation wave. The velocity span of the residual 
profile variation is determined by a convolution of the intrinsic stel- 
lar line profile with the pulsation velocity field. If the former is large 
compared to the pulsation velocity amplitude, as is the case for REE 
lines in roAp stars, the velocity span of LPV is not representative of 
the true pulsation wave speed. This situation is analogous to the be- 
haviour of rapidly rotating non-radial pulsators, where profile vari- 
ation with V p <C v e sin i is observed in the who le velocity range 
from — v c sin i to +v c sin i dSchriivers et al .119971) . Thus, we assert 
that observations of 7 Equ and other roAp stars give no evidence 
of supersonic pu l sation velocity speeds and therefore the model of 
Shibahashi et alj d2004b cannot be correct. 

However, it is still instructive to investigate the consequences 
of a non-sinusoidal velocity variation for oblique pu lsators. In do- 
ing this, we have followed Shibahashi eta[| d2004l) and replaced 
the exp (iwi) term of Eq. by n -1 exp (inivt). Addi- 

tional moderate smoothing along the time axis is required to get 
rid of wiggles introduced by a finite number of terms in the lat- 
ter expression. The line profile and moment variability resulting 
from a saw-tooth time dependence of the pulsation velocity is pre- 
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sented in Fig. [16] Evidently, the residual LPV shows a blue-to- 
red shift, which is superficially similar to the observed profile be- 
haviour. However, the time series of residuals also shows a promi- 
nent gap due to the very fast evolution of the profile shape at phase 
0. This is essentially the signature of the shock wave postulated 
by Shibaha shi et alj d2004l) . No such rapid changes of line profiles 
have been reported for any roAp star. Also, in contradiction with 
observations, the variation of the line profile moments predicted by 
the shock wave model contains a large harmonic contribution. The 
observed phase relation of RV and line width is not reproduced 
either: the maximum line width coincides with the velocity ex- 
tremum, whereas observations point to a 0.25 phase offset between 
the two curves. 

In summary, our theoretical calculations suggest that, even 
putting aside the dubious foundation of the shock wave model, the 
line profile and moment variability predicted by this theory dis- 
agrees with ob servations. In our opinio n, this proves that the model 
formulated by Shiba hashi et al.1 d2004l) does not provide a viable 
solution to the puzzle of LPV in roAp stars. 
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Figure 17. Comparison of very high resolution observation of 7 Equ (thin 
curve) with magnetic spectrum synthesis calculations. Zeeman resolved 
profiles of the iron-peak lines (right panel) are well reproduced with a 
macroturbulent broadening of 2 kms -1 (thick solid curve). In contrast, a 
satisfactory fit to the broad wings of the Ndm 6145 A line (left panel) is 
only possible for Vmacro 10 kms -1 (dashed curve). 



6 DISCUSSION 

Our survey of pulsational line profile variability has, for the 
first time, characterized in detail the rapid fluctuations of REE 
lines in a representative sample of slowly rotating roAp stars. 
For a subset of lines in each star studied, we demonstrated 
the presence of asymmetri c LPV of the type first reported by 
iKochukhov & Rvabchikovj ( feOOlal) for 7 Equ. We point out the 
close relation of this variability to the high-amplitude oscillation 
of the second moment. Both the asymmetric, blue-to-red shifting 
features in the residual profile time series, and the single-wave 
line width modulation with pulsation phase, increase rapidly with 
height and reach a maximum in the uppermost atmospheric layers, 
probed by strong doubly ionized REE lines. We show that this type 
of LPV is incompatible with the standard oblique pulsator model, 
which attributes pulsational spectroscopic variability to the veloc- 
ity perturbations described by the I = 1, m = spherical har- 
monic. Extensive line profile calculations reveal that none of the 
recent modifications and extension s of the oblique pulsator frame- 



work ( iBigot & Dziembowskill2002l ; ISaio & Gautschvll2004l : ICunhal 



2006), nor the exotic shock wave model o f lShibahashi et alj 120041) . 
are capable of explaining the LPV observed in sharp-lined roAp 
stars. Neveretheless, we show that theoretical line profile variation 
can be brought into remarkable agreement with experimental data 
if a pulsational modulation of the width of the local intrinsic spec- 
tral line profile is incorporated into the model. This new oblique 
pulsator framework is based on dipolar modes, and thus does not 
contradict observational evidence of the presence of this type of 
non-radial oscillations, accumulated by numerous previous spec- 
troscopic and photometric studies of roAp stars. At the same time, 
a combination of the velocity and line width variability solves the 
puzzle of asymmetric LPV, explains why this line profile behaviour 
is clearly detected only in sharp-lined roAp stars, and provides a 
satisfactory interpretation of the line profile moments and the mor- 
phology of the standard deviation spectrum. The empirical success 
of our line width modulation hypothesis notwithstanding, the im- 
portant question of the physical mechanism responsible for the pul- 
sational alternation of the width of spectral lines has not yet been 
addressed. Here we attempt to shed light on this problem. 

High-resolution, high-S/N observations of cool Ap stars re- 
veal a number of line profile and line strengths anomalies, usu- 



ally attributed to the combined effects of chemical stratifica- 
tion, a strong magnetic fi eld, and non-solar chemical composition 
dRvabchikova et al.|[2002h . However, the significant broadening of 
strong doubly ionized REE lines in slowly rotating cool Ap stars 
has received no satisfactory explanation so far. This problem is il- 
lustrated in Fig. [T7] for the Ndm 6145 A line in 7 Equ. We com- 
pare the highest resolution observation of this star available in the 
ESQ Archive (coadd ed A/AA = 220,000 spectra obtained by 
ISavanov et al. I ( 120051) ) with the magnetic spectrum synthesis cal- 
culations based on the stell ar parameters, chemical abundances 
and magnetic field model o f iRvabchikova et al.1 (2002). The addi- 
tional broadening required to fit the Zeeman resolved components 
of the Fe II and Cr II lines is approximately 2 km s -1 . This is much 
smaller than V maC ro ^ 10 kms -1 required to roughly reproduce 
the profile of the Nd III line. Such a high Gaussian smearing washes 
out the Zeeman structure in the Fe-peak lines and, therefore, con- 
tradicts observations of spectral lines formed in the lower atmo- 
sphere of 7 Equ. The factor of five discrepancy in the best-fitting 
macroturbulence derived for iron-peak elements and Nd III cannot 
be ascribed to isotope sp litting or hyperfin e structure effects (those 
are small for Nd III, see iDolk et all d2002l) ) and seems to reflect a 
genuine increase of the isotropic velocity component towards the 
high-lying REE cloud. 

Unusually strong broadening of pulsating REE lines w as first 
discussed for 7 Equ by Kochukhov & Rvabchikov^ d2001al) . In that 
paper we advocated a high-^, non-axisymmetric mode in 7 Equ and 
attributed the large broadening of Pr III and Nd III lines to the unre- 
solved, large-amplitude pulsational velocity field. Our new inves- 
tigation of this and other roAp stars suggests a different picture: 
pulsation modes are axisymmetric with 1=1, and line width mod- 
ulation is invoked to explain the observed LPV. This model implies 
local pulsational perturbations with amplitudes not exceeding 1- 
2 kms -1 for the REE line-forming layer, which leaves no room 
for ~ 10 km s _1 pulsational broadening in the time-averaged spec- 
trum. Comparison of the pulsational amplitude and line broadening 
for different time-resolved spectroscopic data sets of 7 Equ and for 
different roAp stars also argues against a direct relation between 
the average line width and non-radial pulsation. We also find no 
differences between the REE line broadening for time-averaged 
spectra obtained from time series of 7 Equ characterized by an 
overall RV amplitude that is different by a factor of five. Many 
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roAp stars have relatively low (<100ms ) pulsation amplitudes 
(HD 137949, HD 166473, HD 176232) but, nevertheless, show very 
broad REE lines. Finally, apart from the small line width modula- 
tion discussed in this paper, no systematic reduction of the REE 
line widths is observed in time-resolved spectra. 

Thus, effects other than pulsational variability must be re- 
sponsible for the inexplicable broadening of REE lines. Turbulence 
caused by convective instability is an obvious candidate for expla- 
nations of roAp-star observations. The location of cool Ap stars in 
the H-R diagramm coincides with Am stars, for which large broad- 
ening of strong line s is observed and interpreted as a signature of 
convective motions dLandstreetl 1 1 99 a ) . On the other hand, the in- 
tense magnetic fields of Ap stars are likely to modify and, possi- 
bly, suppress H and He convection zones. Ubiquitous signatures 
of chemical stratification also indicate a lack of convective mix- 
ing in the lower atmosphere of cool Ap stars, though some mixing 
is required to reduce the di screpancy between the oretical diffusion 
models and observations dLeBlanc & Monir]|2004h . The message 
of these modern theoretical and observational results is that mag- 
netic field-convection interaction in cool Ap stars is far more com- 
plex than previously throught and that a common assumption of the 
full suppression of convection is likely to be an oversimplification. 

With this background, turbulence related to convective or 
pulsationally-driven instability is the most plausible explanation of 
the anomalous broadening of REE lines. We speculate that some 
cool Ap stars, including the roAp pulsators studied in our pa- 
per, possess a turbulence zone in the upper atmospheric layers, 
roughly at the height of the REE-rich cloud. Large vertical gra- 
dients of chemical composition in the intermediate atmosphere 
(—3.5 Si logT5ooo ^ —0.5) exclude significant mixing. Therefore, 
turbulence should appear only in layers above log T5000 ~ — 3.5 
to —4.0. We find that strong REE lines often display excessive 
width of the outer wings, which is difficult to reproduce assum- 
ing a height-independent broadening parameter. Possibly, line cores 
sample layers with a weaker turbulence - a situation consistent with 
a decrease of convection efficiency with height. 

At this stage we can offer no complete physical theory ex- 
plaining the appearance of the turbulence zone in the upper atmo- 
spheres of cool Ap stars. One reason for the convective instabil- 
ity to develop may be a non-standard temperature structure of Ap- 
star atmosp heres. Empirica l evide nce for this kind of anomaly was 
found bv lKochukhov et alj |2002). By fitting the peculiar shape of 
the core- wing transition of the hydrogen Balmer lines, they derived 
semi-empirical models, characterized by a nearly flat temperature 
profile above logTfcooo ~ — 1.0 followed by a steep drop of tem- 
perature at logT5ooo ~ —4.0. This large temperature gradient may 
lead to convective instability very close to the inferred location of 
the turbulence zone needed to explain broadening of REE lines. 

What could be the response of the hypothetical turbulent layer 
to the non-radial p-mode oscillations? We know that the amplitude 
of conve ctive motions grow s as the mean atmospheric pressure in- 
creases faeiter et al.l|20 02). Pulsational modulation of convection 
efficiency was invo ked to explain observations for several classes 
of pu lsating stars jFokin. Gillet & Breitfellned 1 19961 : iGillet et all 
2002). One usually infers that convective motions become stronger 
during the compression phase of the oscillation cycle. Conversly, 
convection is less intense during the expansion of the pulsating star. 
This is precisely the relationship between the line width and radius 
variation that we have deduced for roAp stars. Consequently, it is 
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reasonable to suggest that pulsations modulate the efficiency of tur- 
bulence in the upper atmospheres of roAp stars, and this leads to the 
observed line width and line profile variation. 

We hope that our interpretation of the pulsational LPV in roAp 
stars and the possible discovery of a new turbulence zone will stim- 
ulate new detailed magnetohydrodynamic modelling of the inter- 
action between convection, pulsation and magnetic field in cool 
Ap stars. Three-dimension al radiative hydrodynam ic simulations 
of A-star convection (e.g.. iKochukhov et alj|2006h should be ex- 
tended to include large-scale fossil magnetic fields. Variable lower 
boundary condition can be used to simulate the effects of high- 
overtone p-modes. On the other hand, a NLTE semi-empirical anal- 
ysis of hydrogen lines, coupled with model atmosphere calcula- 
tions takin g into account chemical stratification and the magn etic 
field JShulvak et ai1l2004l : IKochukhov. Khan & Shulvakll2005h. is 
necess ary to verify and extended the results of Kochukho v et al.l 
d2002l) . 
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